I am so excited to get an opportunity to be into R. I am expecting to learn how to hand data and visulaze the results of analysis through R. One of my colleagues had let me know about this course.
The URL for my git repository as follow:
https://github.com/Donggyu-Kam/IODS-project
The data is the student’s attitude toward statistics.The data consists of information about students’ ages, students’ gender, three approaches, students’ attitude toward statistics and exam points. The oberservations whare the exam points are zeor were removed from data. And thus, the data has 166 observations and 7 variables. You can figure out three variables to affect the exam points through linear regression models and finally you will figure out the most significant variable to exampoints
#Before handling data, you should install two packages:“ggplot2” and“GGally”.
By drawing plots illustrating the relationships between the variables and the normality of the erros distributions, it will find out what variables are closely related to exam points.
According to the plot “p”, three variables “attitude”, “stra” and “surf” were chosen due to the strongest correlation value as explanatory variables.
After fitting each variable to a linear regression modle, “attitude” shows postive relationship but the others show week or negative relationships. So, the linear regression model is re-fitted with three variables together.
In this model, “Stra” and “Surf” are not not significantly related to the exam points. So, the model is re-fitted again witout the two variables
Fit the model again with the “attitude”
By analyzing the residuals and errors of the model, we can test the model validation. The variance of errors in residulas over fitted values looks constant without any patterns but there are 4 outliers indicating that the model only with “attitude” cannot explain the 4 ouliters of exam points. In the normality of errors in QQ plot, the errors are normaly distributed. In leverage of observations, standardized residuals are scattered around 0 and there seems no impact of outliers
The data is about student acheivement in secondary education of two Portugues schools. Two datasets are available for the performance in two subjects: Math and Portuguese language.
P. Cortez and A. Silva. Using Data Mining to Predict Secondary School Student Performance. In A. Brito and J. Teixeira Eds., Proceedings of 5th FUture BUsiness TEChnology Conference (FUBUTEC 2008) pp. 5-12, Porto, Portugal, April, 2008, EUROSIS, ISBN 978-9077381-39-7.[link]http://www3.dsi.uminho.pt/pcortez/student.pdf
the data is about student acheivement in secondary education of two Portugues schools according to the student grades, demographic and school-related features. Two datasets are available for the performance in two subjects: Math and Portuguese language. the datasets consists of 35 variables and 382 observations.I firstly beging with installing some packages helps analysis. And thus, we can take a glimpse at the datasets and what variables are available through column names on the datasets.
library(dplyr);library(ggplot2);library(tidyr);library(boot)
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
alc<-read.csv("Z:\\IODS-project\\data\\alc.csv")
glimpse(alc)
## Observations: 382
## Variables: 36
## $ X <int> 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, ...
## $ school <fct> GP, GP, GP, GP, GP, GP, GP, GP, GP, GP, GP, GP, GP,...
## $ sex <fct> F, F, F, F, F, M, M, F, M, M, F, F, M, M, M, F, F, ...
## $ age <int> 18, 17, 15, 15, 16, 16, 16, 17, 15, 15, 15, 15, 15,...
## $ address <fct> U, U, U, U, U, U, U, U, U, U, U, U, U, U, U, U, U, ...
## $ famsize <fct> GT3, GT3, LE3, GT3, GT3, LE3, LE3, GT3, LE3, GT3, G...
## $ Pstatus <fct> A, T, T, T, T, T, T, A, A, T, T, T, T, T, A, T, T, ...
## $ Medu <int> 4, 1, 1, 4, 3, 4, 2, 4, 3, 3, 4, 2, 4, 4, 2, 4, 4, ...
## $ Fedu <int> 4, 1, 1, 2, 3, 3, 2, 4, 2, 4, 4, 1, 4, 3, 2, 4, 4, ...
## $ Mjob <fct> at_home, at_home, at_home, health, other, services,...
## $ Fjob <fct> teacher, other, other, services, other, other, othe...
## $ reason <fct> course, course, other, home, home, reputation, home...
## $ nursery <fct> yes, no, yes, yes, yes, yes, yes, yes, yes, yes, ye...
## $ internet <fct> no, yes, yes, yes, no, yes, yes, no, yes, yes, yes,...
## $ guardian <fct> mother, father, mother, mother, father, mother, mot...
## $ traveltime <int> 2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 1, 1, ...
## $ studytime <int> 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 3, 1, 2, 3, 1, 3, ...
## $ failures <int> 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...
## $ schoolsup <fct> yes, no, yes, no, no, no, no, yes, no, no, no, no, ...
## $ famsup <fct> no, yes, no, yes, yes, yes, no, yes, yes, yes, yes,...
## $ paid <fct> no, no, yes, yes, yes, yes, no, no, yes, yes, yes, ...
## $ activities <fct> no, no, no, yes, no, yes, no, no, no, yes, no, yes,...
## $ higher <fct> yes, yes, yes, yes, yes, yes, yes, yes, yes, yes, y...
## $ romantic <fct> no, no, no, yes, no, no, no, no, no, no, no, no, no...
## $ famrel <int> 4, 5, 4, 3, 4, 5, 4, 4, 4, 5, 3, 5, 4, 5, 4, 4, 3, ...
## $ freetime <int> 3, 3, 3, 2, 3, 4, 4, 1, 2, 5, 3, 2, 3, 4, 5, 4, 2, ...
## $ goout <int> 4, 3, 2, 2, 2, 2, 4, 4, 2, 1, 3, 2, 3, 3, 2, 4, 3, ...
## $ Dalc <int> 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...
## $ Walc <int> 1, 1, 3, 1, 2, 2, 1, 1, 1, 1, 2, 1, 3, 2, 1, 2, 2, ...
## $ health <int> 3, 3, 3, 5, 5, 5, 3, 1, 1, 5, 2, 4, 5, 3, 3, 2, 2, ...
## $ absences <int> 6, 4, 10, 2, 4, 10, 0, 6, 0, 0, 0, 4, 2, 2, 0, 4, 6...
## $ G1 <int> 5, 5, 7, 15, 6, 15, 12, 6, 16, 14, 10, 10, 14, 10, ...
## $ G2 <int> 6, 5, 8, 14, 10, 15, 12, 5, 18, 15, 8, 12, 14, 10, ...
## $ G3 <int> 6, 6, 10, 15, 10, 15, 11, 6, 19, 15, 9, 12, 14, 11,...
## $ alc_use <dbl> 1.0, 1.0, 2.5, 1.0, 1.5, 1.5, 1.0, 1.0, 1.0, 1.0, 1...
## $ high_use <lgl> TRUE, TRUE, FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, TR...
gather(alc)%>% glimpse
## Warning: attributes are not identical across measure variables;
## they will be dropped
## Observations: 13,752
## Variables: 2
## $ key <chr> "X", "X", "X", "X", "X", "X", "X", "X", "X", "X", "X", "...
## $ value <chr> "1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11",...
gather(alc)%>%ggplot(aes(value))+facet_wrap("key", scales="free")+geom_bar()
## Warning: attributes are not identical across measure variables;
## they will be dropped
In order to study the relationship between high and low alcohol consumptions and variables, i am going to work folloing to hypothesis;
I will take a look at each relationship followring to the hypothesis one by one p<-ggpairs(alc, mapping = aes(col=sex, alpha = 0.3), lower = list(combo = wrap(“facethist”, bins = 20))) p
g1 <- ggplot(data = alc, aes(x = alc_use, fill = sex))
g1 + geom_bar()+facet_wrap("sex") + ggtitle("alohol consumption by sex")
g11<-ggplot(data=alc, aes(x=high_use, fill=sex))
g11+geom_bar()+facet_wrap("sex") + ggtitle("high and low alcohol consumption by sex")
g2 <- ggplot(data=alc, aes(x=high_use, y= G1, col=sex))
g2+geom_boxplot()+ggtitle("G1 by alcohol consumption and sex")
g3<-ggplot(alc, aes(x=high_use, y=G2, col=sex))
g3+geom_boxplot()+ggtitle("G2 by acohol consumption and sex")
g4<-ggplot(alc, aes(x=high_use, y=G3, col=sex))
g4+geom_boxplot()+ggtitle("G3 by alcohol consumption and sex")
As a result of exploration of relationship with 4 variable to the high and low alcohol relationship, it seems that femailes consume more alcohol than males do. However, it is presumed that the more males consum alcohol, the less grades males achieve.
From now, in order to interpret relationship between variables presumably relating to alcohol consumptions and alcohol consumption, we will deal with logistaic regression.
m<-glm(high_use~sex+G1+G2+G3, data=alc, family="binomial")
summary(m)
##
## Call:
## glm(formula = high_use ~ sex + G1 + G2 + G3, family = "binomial",
## data = alc)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.7050 -1.1393 0.8236 1.0275 1.4369
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.124551 0.377648 -0.330 0.742
## sexM -0.929752 0.215442 -4.316 1.59e-05 ***
## G1 0.100421 0.063556 1.580 0.114
## G2 0.005342 0.076244 0.070 0.944
## G3 -0.032884 0.053600 -0.614 0.540
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 524.94 on 381 degrees of freedom
## Residual deviance: 502.43 on 377 degrees of freedom
## AIC: 512.43
##
## Number of Fisher Scoring iterations: 4
coef(m)
## (Intercept) sexM G1 G2 G3
## -0.124550963 -0.929752246 0.100420779 0.005341903 -0.032884237
OR<-coef(m)%>%exp
OR
## (Intercept) sexM G1 G2 G3
## 0.8828933 0.3946515 1.1056360 1.0053562 0.9676506
CI<-exp(confint(m))
## Waiting for profiling to be done...
CI
## 2.5 % 97.5 %
## (Intercept) 0.4202206 1.8520181
## sexM 0.2575779 0.5999088
## G1 0.9774251 1.2553494
## G2 0.8649561 1.1681070
## G3 0.8686826 1.0737634
cbind(OR, CI)
## OR 2.5 % 97.5 %
## (Intercept) 0.8828933 0.4202206 1.8520181
## sexM 0.3946515 0.2575779 0.5999088
## G1 1.1056360 0.9774251 1.2553494
## G2 1.0053562 0.8649561 1.1681070
## G3 0.9676506 0.8686826 1.0737634
The model indicates taht males are more significantly related to high alcohol consumptions.However, the other 3 variables are not significantly related to high alcohol consumptions. In odds rates, the odds rates for G2 and G3 are more than 1.
probabilities<-predict(m, type="response")
alc<-mutate(alc,probability=probabilities)
alc<-mutate(alc,prediction=probability>0.5)
select(alc,G2,G3, high_use, probability, prediction)%>% tail(10)
## G2 G3 high_use probability prediction
## 373 5 0 FALSE 0.3953044 FALSE
## 374 5 5 FALSE 0.3567518 FALSE
## 375 9 10 FALSE 0.6454027 TRUE
## 376 5 6 FALSE 0.5762451 TRUE
## 377 5 0 TRUE 0.6468232 TRUE
## 378 9 8 TRUE 0.5898627 TRUE
## 379 5 0 TRUE 0.6235593 TRUE
## 380 5 0 TRUE 0.6235593 TRUE
## 381 16 16 FALSE 0.4777423 FALSE
## 382 12 10 FALSE 0.4466009 FALSE
g<-ggplot(alc,(aes(x=probability, y=high_use, col=prediction)))
g+geom_point()
table(high_use=alc$high_use, prediction=alc$prediction)%>%prop.table()%>%addmargins()
## prediction
## high_use FALSE TRUE Sum
## FALSE 0.2225131 0.2225131 0.4450262
## TRUE 0.1492147 0.4057592 0.5549738
## Sum 0.3717277 0.6282723 1.0000000
loss_func<-function(class, prob){
n_wrong <-abs(class - prob) > 0.5
mean(n_wrong)
}
loss_func(class = alc$high_use, prob = alc$probability)
## [1] 0.3717277
Through exploring the training error, the incorrection of prediction is almost 50%.
cv<-cv.glm(data=alc, cost=loss_func, glmfit = m, K=10)
cv$delta[1]
## [1] 0.382199
library(MASS);library(tidyverse);library(corrplot)
##
## Attaching package: 'MASS'
## The following object is masked from 'package:dplyr':
##
## select
## -- Attaching packages -------------------------------------------------------------------------------------- tidyverse 1.2.1 --
## v tibble 2.1.3 v purrr 0.3.3
## v readr 1.3.1 v stringr 1.4.0
## v tibble 2.1.3 v forcats 0.4.0
## -- Conflicts ----------------------------------------------------------------------------------------- tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag() masks stats::lag()
## x MASS::select() masks dplyr::select()
## corrplot 0.84 loaded
data("Boston")
The data has 506 observations and 14 variables. Data is about housing values in suburbs of Boston.
str(Boston)
## 'data.frame': 506 obs. of 14 variables:
## $ crim : num 0.00632 0.02731 0.02729 0.03237 0.06905 ...
## $ zn : num 18 0 0 0 0 0 12.5 12.5 12.5 12.5 ...
## $ indus : num 2.31 7.07 7.07 2.18 2.18 2.18 7.87 7.87 7.87 7.87 ...
## $ chas : int 0 0 0 0 0 0 0 0 0 0 ...
## $ nox : num 0.538 0.469 0.469 0.458 0.458 0.458 0.524 0.524 0.524 0.524 ...
## $ rm : num 6.58 6.42 7.18 7 7.15 ...
## $ age : num 65.2 78.9 61.1 45.8 54.2 58.7 66.6 96.1 100 85.9 ...
## $ dis : num 4.09 4.97 4.97 6.06 6.06 ...
## $ rad : int 1 2 2 3 3 3 5 5 5 5 ...
## $ tax : num 296 242 242 222 222 222 311 311 311 311 ...
## $ ptratio: num 15.3 17.8 17.8 18.7 18.7 18.7 15.2 15.2 15.2 15.2 ...
## $ black : num 397 397 393 395 397 ...
## $ lstat : num 4.98 9.14 4.03 2.94 5.33 ...
## $ medv : num 24 21.6 34.7 33.4 36.2 28.7 22.9 27.1 16.5 18.9 ...
summary(Boston)
## crim zn indus chas
## Min. : 0.00632 Min. : 0.00 Min. : 0.46 Min. :0.00000
## 1st Qu.: 0.08204 1st Qu.: 0.00 1st Qu.: 5.19 1st Qu.:0.00000
## Median : 0.25651 Median : 0.00 Median : 9.69 Median :0.00000
## Mean : 3.61352 Mean : 11.36 Mean :11.14 Mean :0.06917
## 3rd Qu.: 3.67708 3rd Qu.: 12.50 3rd Qu.:18.10 3rd Qu.:0.00000
## Max. :88.97620 Max. :100.00 Max. :27.74 Max. :1.00000
## nox rm age dis
## Min. :0.3850 Min. :3.561 Min. : 2.90 Min. : 1.130
## 1st Qu.:0.4490 1st Qu.:5.886 1st Qu.: 45.02 1st Qu.: 2.100
## Median :0.5380 Median :6.208 Median : 77.50 Median : 3.207
## Mean :0.5547 Mean :6.285 Mean : 68.57 Mean : 3.795
## 3rd Qu.:0.6240 3rd Qu.:6.623 3rd Qu.: 94.08 3rd Qu.: 5.188
## Max. :0.8710 Max. :8.780 Max. :100.00 Max. :12.127
## rad tax ptratio black
## Min. : 1.000 Min. :187.0 Min. :12.60 Min. : 0.32
## 1st Qu.: 4.000 1st Qu.:279.0 1st Qu.:17.40 1st Qu.:375.38
## Median : 5.000 Median :330.0 Median :19.05 Median :391.44
## Mean : 9.549 Mean :408.2 Mean :18.46 Mean :356.67
## 3rd Qu.:24.000 3rd Qu.:666.0 3rd Qu.:20.20 3rd Qu.:396.23
## Max. :24.000 Max. :711.0 Max. :22.00 Max. :396.90
## lstat medv
## Min. : 1.73 Min. : 5.00
## 1st Qu.: 6.95 1st Qu.:17.02
## Median :11.36 Median :21.20
## Mean :12.65 Mean :22.53
## 3rd Qu.:16.95 3rd Qu.:25.00
## Max. :37.97 Max. :50.00
pairs(Boston)
cor_matrix<-cor(Boston)
cor_matrix%>%round(2)
## crim zn indus chas nox rm age dis rad tax
## crim 1.00 -0.20 0.41 -0.06 0.42 -0.22 0.35 -0.38 0.63 0.58
## zn -0.20 1.00 -0.53 -0.04 -0.52 0.31 -0.57 0.66 -0.31 -0.31
## indus 0.41 -0.53 1.00 0.06 0.76 -0.39 0.64 -0.71 0.60 0.72
## chas -0.06 -0.04 0.06 1.00 0.09 0.09 0.09 -0.10 -0.01 -0.04
## nox 0.42 -0.52 0.76 0.09 1.00 -0.30 0.73 -0.77 0.61 0.67
## rm -0.22 0.31 -0.39 0.09 -0.30 1.00 -0.24 0.21 -0.21 -0.29
## age 0.35 -0.57 0.64 0.09 0.73 -0.24 1.00 -0.75 0.46 0.51
## dis -0.38 0.66 -0.71 -0.10 -0.77 0.21 -0.75 1.00 -0.49 -0.53
## rad 0.63 -0.31 0.60 -0.01 0.61 -0.21 0.46 -0.49 1.00 0.91
## tax 0.58 -0.31 0.72 -0.04 0.67 -0.29 0.51 -0.53 0.91 1.00
## ptratio 0.29 -0.39 0.38 -0.12 0.19 -0.36 0.26 -0.23 0.46 0.46
## black -0.39 0.18 -0.36 0.05 -0.38 0.13 -0.27 0.29 -0.44 -0.44
## lstat 0.46 -0.41 0.60 -0.05 0.59 -0.61 0.60 -0.50 0.49 0.54
## medv -0.39 0.36 -0.48 0.18 -0.43 0.70 -0.38 0.25 -0.38 -0.47
## ptratio black lstat medv
## crim 0.29 -0.39 0.46 -0.39
## zn -0.39 0.18 -0.41 0.36
## indus 0.38 -0.36 0.60 -0.48
## chas -0.12 0.05 -0.05 0.18
## nox 0.19 -0.38 0.59 -0.43
## rm -0.36 0.13 -0.61 0.70
## age 0.26 -0.27 0.60 -0.38
## dis -0.23 0.29 -0.50 0.25
## rad 0.46 -0.44 0.49 -0.38
## tax 0.46 -0.44 0.54 -0.47
## ptratio 1.00 -0.18 0.37 -0.51
## black -0.18 1.00 -0.37 0.33
## lstat 0.37 -0.37 1.00 -0.74
## medv -0.51 0.33 -0.74 1.00
corrplot(cor_matrix, method="circle", type="upper", cl.pos="b", tl.pos="d", tl.cex=0.6)
There is strongly positive correlation bewtween “rad” and “tax” and stongly negative correlation between “age” and “dis” and between “lstat” and “medv”.
boston_scaled<-scale(Boston)
summary(boston_scaled)
## crim zn indus
## Min. :-0.419367 Min. :-0.48724 Min. :-1.5563
## 1st Qu.:-0.410563 1st Qu.:-0.48724 1st Qu.:-0.8668
## Median :-0.390280 Median :-0.48724 Median :-0.2109
## Mean : 0.000000 Mean : 0.00000 Mean : 0.0000
## 3rd Qu.: 0.007389 3rd Qu.: 0.04872 3rd Qu.: 1.0150
## Max. : 9.924110 Max. : 3.80047 Max. : 2.4202
## chas nox rm age
## Min. :-0.2723 Min. :-1.4644 Min. :-3.8764 Min. :-2.3331
## 1st Qu.:-0.2723 1st Qu.:-0.9121 1st Qu.:-0.5681 1st Qu.:-0.8366
## Median :-0.2723 Median :-0.1441 Median :-0.1084 Median : 0.3171
## Mean : 0.0000 Mean : 0.0000 Mean : 0.0000 Mean : 0.0000
## 3rd Qu.:-0.2723 3rd Qu.: 0.5981 3rd Qu.: 0.4823 3rd Qu.: 0.9059
## Max. : 3.6648 Max. : 2.7296 Max. : 3.5515 Max. : 1.1164
## dis rad tax ptratio
## Min. :-1.2658 Min. :-0.9819 Min. :-1.3127 Min. :-2.7047
## 1st Qu.:-0.8049 1st Qu.:-0.6373 1st Qu.:-0.7668 1st Qu.:-0.4876
## Median :-0.2790 Median :-0.5225 Median :-0.4642 Median : 0.2746
## Mean : 0.0000 Mean : 0.0000 Mean : 0.0000 Mean : 0.0000
## 3rd Qu.: 0.6617 3rd Qu.: 1.6596 3rd Qu.: 1.5294 3rd Qu.: 0.8058
## Max. : 3.9566 Max. : 1.6596 Max. : 1.7964 Max. : 1.6372
## black lstat medv
## Min. :-3.9033 Min. :-1.5296 Min. :-1.9063
## 1st Qu.: 0.2049 1st Qu.:-0.7986 1st Qu.:-0.5989
## Median : 0.3808 Median :-0.1811 Median :-0.1449
## Mean : 0.0000 Mean : 0.0000 Mean : 0.0000
## 3rd Qu.: 0.4332 3rd Qu.: 0.6024 3rd Qu.: 0.2683
## Max. : 0.4406 Max. : 3.5453 Max. : 2.9865
class(boston_scaled)
## [1] "matrix"
boston_scaled<-as.data.frame(boston_scaled)
summary("crim")
## Length Class Mode
## 1 character character
bins<-quantile(boston_scaled$crim)
bins
## 0% 25% 50% 75% 100%
## -0.419366929 -0.410563278 -0.390280295 0.007389247 9.924109610
crime<-cut(boston_scaled$crim, breaks=bins, include.lowest=TRUE, label=c("low","med_low","med_high", "high"))
crime
## [1] low low low low low low med_low
## [8] med_low med_low med_low med_low med_low med_low med_high
## [15] med_high med_high med_high med_high med_high med_high med_high
## [22] med_high med_high med_high med_high med_high med_high med_high
## [29] med_high med_high med_high med_high med_high med_high med_high
## [36] low med_low low med_low low low med_low
## [43] med_low med_low med_low med_low med_low med_low med_low
## [50] med_low med_low low low low low low
## [57] low low med_low med_low med_low med_low med_low
## [64] med_low low low low low med_low med_low
## [71] med_low med_low med_low med_low low med_low med_low
## [78] med_low low med_low low low low low
## [85] low low low low low low low
## [92] low low low low med_low med_low med_low
## [99] low low med_low med_low med_low med_low med_low
## [106] med_low med_low med_low med_low med_high med_low med_low
## [113] med_low med_low med_low med_low med_low med_low med_low
## [120] med_low low low med_low med_low med_low med_low
## [127] med_high med_high med_high med_high med_high med_high med_high
## [134] med_high med_high med_high med_high med_high med_low med_high
## [141] med_high med_high med_high high med_high med_high med_high
## [148] med_high med_high med_high med_high med_high med_high med_high
## [155] med_high med_high med_high med_high med_high med_high med_high
## [162] med_high med_high med_high med_high med_high med_high med_high
## [169] med_high med_high med_high med_high med_low med_low med_low
## [176] low low low low low low low
## [183] med_low med_low med_low low low low med_low
## [190] med_low med_low low med_low low low low
## [197] low low low low low low low
## [204] low low med_low med_low med_low med_low med_high
## [211] med_low med_high med_low med_low med_high med_low low
## [218] low med_low med_low med_high med_high med_high med_high
## [225] med_high med_high med_high med_high med_high med_high med_high
## [232] med_high med_high med_high med_high med_high med_high med_high
## [239] med_low med_low med_low med_low med_low med_low med_low
## [246] med_low med_high med_low med_low med_low med_low med_low
## [253] med_low med_high low low low med_high med_high
## [260] med_high med_high med_high med_high med_high med_high med_high
## [267] med_high med_high med_high med_low med_high med_low med_low
## [274] med_low low med_low med_low low low med_low
## [281] low low low low low low low
## [288] low low low low low low med_low
## [295] low med_low low med_low low low low
## [302] low med_low med_low low low low low
## [309] med_high med_high med_high med_high med_high med_high med_high
## [316] med_low med_high med_low med_high med_high med_low med_low
## [323] med_high med_high med_high med_low med_high med_low low
## [330] low low low low low low low
## [337] low low low low low low low
## [344] low low low low low low low
## [351] low low low low low med_low high
## [358] high high high high high high high
## [365] med_high high high high high high high
## [372] high high high high high high high
## [379] high high high high high high high
## [386] high high high high high high high
## [393] high high high high high high high
## [400] high high high high high high high
## [407] high high high high high high high
## [414] high high high high high high high
## [421] high high high high high high high
## [428] high high high high high high high
## [435] high high high high high high high
## [442] high high high high high high high
## [449] high high high high high high high
## [456] high high high high high high high
## [463] high high high med_high high high high
## [470] high high high med_high high high high
## [477] high high high high high high high
## [484] med_high med_high med_high high high med_low med_low
## [491] med_low med_low med_low med_low med_high med_low med_high
## [498] med_high med_low med_low med_low low low low
## [505] med_low low
## Levels: low med_low med_high high
table(crime)
## crime
## low med_low med_high high
## 127 126 126 127
boston_scaled<-dplyr::select(boston_scaled,-crim)
boston_scaled
## zn indus chas nox rm
## 1 0.28454827 -1.28663623 -0.2723291 -0.14407485 0.413262920
## 2 -0.48724019 -0.59279438 -0.2723291 -0.73953036 0.194082387
## 3 -0.48724019 -0.59279438 -0.2723291 -0.73953036 1.281445551
## 4 -0.48724019 -1.30558569 -0.2723291 -0.83445805 1.015297761
## 5 -0.48724019 -1.30558569 -0.2723291 -0.83445805 1.227362043
## 6 -0.48724019 -1.30558569 -0.2723291 -0.83445805 0.206891639
## 7 0.04872402 -0.47618230 -0.2723291 -0.26489191 -0.388026950
## 8 0.04872402 -0.47618230 -0.2723291 -0.26489191 -0.160306916
## 9 0.04872402 -0.47618230 -0.2723291 -0.26489191 -0.930285282
## 10 0.04872402 -0.47618230 -0.2723291 -0.26489191 -0.399412952
## 11 0.04872402 -0.47618230 -0.2723291 -0.26489191 0.131459378
## 12 0.04872402 -0.47618230 -0.2723291 -0.26489191 -0.392296701
## 13 0.04872402 -0.47618230 -0.2723291 -0.26489191 -0.563086727
## 14 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.477691714
## 15 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.268473932
## 16 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.641365488
## 17 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.497617217
## 18 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.419338455
## 19 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -1.179354069
## 20 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.793653261
## 21 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -1.017103545
## 22 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.454919710
## 23 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.203004422
## 24 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.671253743
## 25 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.513272969
## 26 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.975829289
## 27 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.671253743
## 28 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.338213193
## 29 -0.48724019 -0.43682573 -0.2723291 -0.14407485 0.299402903
## 30 -0.48724019 -0.43682573 -0.2723291 -0.14407485 0.554164692
## 31 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.813578764
## 32 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.302631937
## 33 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.476268464
## 34 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.830657767
## 35 -0.48724019 -0.43682573 -0.2723291 -0.14407485 -0.268473932
## 36 -0.48724019 -0.75459363 -0.2723291 -0.48063666 -0.500463717
## 37 -0.48724019 -0.75459363 -0.2723291 -0.48063666 -0.631402737
## 38 -0.48724019 -0.75459363 -0.2723291 -0.48063666 -0.618593485
## 39 -0.48724019 -0.75459363 -0.2723291 -0.48063666 -0.453496460
## 40 2.72854505 -1.19334657 -0.2723291 -1.09335175 0.441727925
## 41 2.72854505 -1.19334657 -0.2723291 -1.09335175 1.052302267
## 42 -0.48724019 -0.61611679 -0.2723291 -0.92075595 0.690796712
## 43 -0.48724019 -0.61611679 -0.2723291 -0.92075595 -0.164576667
## 44 -0.48724019 -0.61611679 -0.2723291 -0.92075595 -0.104800158
## 45 -0.48724019 -0.61611679 -0.2723291 -0.92075595 -0.306901688
## 46 -0.48724019 -0.61611679 -0.2723291 -0.92075595 -0.857699521
## 47 -0.48724019 -0.61611679 -0.2723291 -0.92075595 -0.709681499
## 48 -0.48724019 -0.61611679 -0.2723291 -0.92075595 -0.362408446
## 49 -0.48724019 -0.61611679 -0.2723291 -0.92075595 -1.260479332
## 50 -0.48724019 -0.61611679 -0.2723291 -0.92075595 -0.971559538
## 51 0.41317968 -0.80123846 -0.2723291 -0.99842406 -0.457766211
## 52 0.41317968 -0.80123846 -0.2723291 -0.99842406 -0.241432178
## 53 0.41317968 -0.80123846 -0.2723291 -0.99842406 0.322174907
## 54 0.41317968 -0.80123846 -0.2723291 -0.99842406 -0.407952453
## 55 2.72854505 -1.04029322 -0.2723291 -1.24868797 -0.564509977
## 56 3.37170210 -1.44552019 -0.2723291 -1.30909650 1.372533565
## 57 3.15731641 -1.51548743 -0.2723291 -1.24868797 0.139998879
## 58 3.80047346 -1.43094368 -0.2723291 -1.24005818 0.756266222
## 59 0.58468822 -0.87557866 -0.2723291 -0.87760700 -0.198734672
## 60 0.58468822 -0.87557866 -0.2723291 -0.87760700 -0.509003218
## 61 0.58468822 -0.87557866 -0.2723291 -0.87760700 -0.773727758
## 62 0.58468822 -0.87557866 -0.2723291 -0.87760700 -0.453496460
## 63 0.58468822 -0.87557866 -0.2723291 -0.87760700 0.243896145
## 64 0.58468822 -0.87557866 -0.2723291 -0.87760700 0.679410711
## 65 0.26310970 -1.42219777 -0.2723291 -1.19604625 1.166162284
## 66 2.94293073 -1.13212523 -0.2723291 -1.35224545 0.007636609
## 67 2.94293073 -1.13212523 -0.2723291 -1.35224545 -0.708258248
## 68 0.04872402 -0.73855947 -0.2723291 -1.25731776 -0.578742479
## 69 0.04872402 -0.73855947 -0.2723291 -1.25731776 -0.982945540
## 70 0.04872402 -0.73855947 -0.2723291 -1.25731776 -0.568779727
## 71 -0.48724019 -0.04763292 -0.2723291 -1.22279860 0.188389387
## 72 -0.48724019 -0.04763292 -0.2723291 -1.22279860 -0.460612711
## 73 -0.48724019 -0.04763292 -0.2723291 -1.22279860 -0.312594689
## 74 -0.48724019 -0.04763292 -0.2723291 -1.22279860 -0.056409650
## 75 -0.48724019 0.24681257 -0.2723291 -1.01568364 -0.016558644
## 76 -0.48724019 0.24681257 -0.2723291 -1.01568364 0.001943608
## 77 -0.48724019 0.24681257 -0.2723291 -1.01568364 -0.008019143
## 78 -0.48724019 0.24681257 -0.2723291 -1.01568364 -0.205850923
## 79 -0.48724019 0.24681257 -0.2723291 -1.01568364 -0.074911903
## 80 -0.48724019 0.24681257 -0.2723291 -1.01568364 -0.584435480
## 81 0.58468822 -0.91493524 -0.2723291 -1.11061133 0.629596953
## 82 0.58468822 -0.91493524 -0.2723291 -1.11061133 0.475885930
## 83 0.58468822 -0.91493524 -0.2723291 -1.11061133 0.024715612
## 84 0.58468822 -0.91493524 -0.2723291 -1.11061133 -0.167423167
## 85 -0.48724019 -0.96886832 -0.2723291 -0.91212616 0.148538381
## 86 -0.48724019 -0.96886832 -0.2723291 -0.91212616 0.491541682
## 87 -0.48724019 -0.96886832 -0.2723291 -0.91212616 -0.383757200
## 88 -0.48724019 -0.96886832 -0.2723291 -0.91212616 -0.232892677
## 89 -0.48724019 -1.12629463 -0.2723291 -0.56693456 1.028107013
## 90 -0.48724019 -1.12629463 -0.2723291 -0.56693456 1.130581029
## 91 -0.48724019 -1.12629463 -0.2723291 -0.56693456 0.188389387
## 92 -0.48724019 -1.12629463 -0.2723291 -0.56693456 0.171310384
## 93 0.71331963 0.56895343 -0.2723291 -0.78267931 0.223970642
## 94 0.71331963 0.56895343 -0.2723291 -0.78267931 -0.104800158
## 95 0.71331963 0.56895343 -0.2723291 -0.78267931 -0.050716649
## 96 -0.48724019 -1.20209248 -0.2723291 -0.94664532 0.484425431
## 97 -0.48724019 -1.20209248 -0.2723291 -0.94664532 -0.173116168
## 98 -0.48724019 -1.20209248 -0.2723291 -0.94664532 2.539598741
## 99 -0.48724019 -1.20209248 -0.2723291 -0.94664532 2.185209437
## 100 -0.48724019 -1.20209248 -0.2723291 -0.94664532 1.610216351
## 101 -0.48724019 -0.37560439 -0.2723291 -0.29941107 0.629596953
## 102 -0.48724019 -0.37560439 -0.2723291 -0.29941107 0.706452465
## 103 -0.48724019 -0.37560439 -0.2723291 -0.29941107 0.171310384
## 104 -0.48724019 -0.37560439 -0.2723291 -0.29941107 -0.210120673
## 105 -0.48724019 -0.37560439 -0.2723291 -0.29941107 -0.167423167
## 106 -0.48724019 -0.37560439 -0.2723291 -0.29941107 -0.617170235
## 107 -0.48724019 -0.37560439 -0.2723291 -0.29941107 -0.638518988
## 108 -0.48724019 -0.37560439 -0.2723291 -0.29941107 -0.224353176
## 109 -0.48724019 -0.37560439 -0.2723291 -0.29941107 0.269514649
## 110 -0.48724019 -0.37560439 -0.2723291 -0.29941107 -0.079181654
## 111 -0.48724019 -0.37560439 -0.2723291 -0.29941107 -0.127572161
## 112 -0.48724019 -0.16424500 -0.2723291 -0.06640675 0.612517950
## 113 -0.48724019 -0.16424500 -0.2723291 -0.06640675 -0.528928721
## 114 -0.48724019 -0.16424500 -0.2723291 -0.06640675 -0.274166933
## 115 -0.48724019 -0.16424500 -0.2723291 -0.06640675 -0.043600398
## 116 -0.48724019 -0.16424500 -0.2723291 -0.06640675 -0.507579968
## 117 -0.48724019 -0.16424500 -0.2723291 -0.06640675 -0.154613915
## 118 -0.48724019 -0.16424500 -0.2723291 -0.06640675 -0.375217698
## 119 -0.48724019 -0.16424500 -0.2723291 -0.06640675 -0.587281980
## 120 -0.48724019 -0.16424500 -0.2723291 -0.06640675 -0.787960260
## 121 -0.48724019 2.11552109 -0.2723291 0.22700611 -0.590128481
## 122 -0.48724019 2.11552109 -0.2723291 0.22700611 -0.399412952
## 123 -0.48724019 2.11552109 -0.2723291 0.22700611 -0.460612711
## 124 -0.48724019 2.11552109 -0.2723291 0.22700611 -0.610053984
## 125 -0.48724019 2.11552109 -0.2723291 0.22700611 -0.577319229
## 126 -0.48724019 2.11552109 -0.2723291 0.22700611 -0.425031456
## 127 -0.48724019 2.11552109 -0.2723291 0.22700611 -0.955903786
## 128 -0.48724019 1.56744433 -0.2723291 0.59808708 -0.842043769
## 129 -0.48724019 1.56744433 -0.2723291 0.59808708 0.208314890
## 130 -0.48724019 1.56744433 -0.2723291 0.59808708 -0.921745781
## 131 -0.48724019 1.56744433 -0.2723291 0.59808708 0.246742645
## 132 -0.48724019 1.56744433 -0.2723291 0.59808708 0.058873617
## 133 -0.48724019 1.56744433 -0.2723291 0.59808708 0.124343127
## 134 -0.48724019 1.56744433 -0.2723291 0.59808708 -0.658444491
## 135 -0.48724019 1.56744433 -0.2723291 0.59808708 -0.750955755
## 136 -0.48724019 1.56744433 -0.2723291 0.59808708 0.071682869
## 137 -0.48724019 1.56744433 -0.2723291 0.59808708 -0.487654465
## 138 -0.48724019 1.56744433 -0.2723291 0.59808708 0.241049645
## 139 -0.48724019 1.56744433 -0.2723291 0.59808708 -0.608630733
## 140 -0.48724019 1.56744433 -0.2723291 0.59808708 -0.190195170
## 141 -0.48724019 1.56744433 -0.2723291 0.59808708 -0.157460416
## 142 -0.48724019 1.56744433 -0.2723291 0.59808708 -1.801314413
## 143 -0.48724019 1.23072696 3.6647712 2.72964520 -1.254786331
## 144 -0.48724019 1.23072696 -0.2723291 2.72964520 -1.162275067
## 145 -0.48724019 1.23072696 -0.2723291 2.72964520 -1.966411438
## 146 -0.48724019 1.23072696 -0.2723291 2.72964520 -0.220083425
## 147 -0.48724019 1.23072696 -0.2723291 2.72964520 -0.934555033
## 148 -0.48724019 1.23072696 -0.2723291 2.72964520 -1.933676683
## 149 -0.48724019 1.23072696 -0.2723291 2.72964520 -1.563631627
## 150 -0.48724019 1.23072696 -0.2723291 2.72964520 -0.978675789
## 151 -0.48724019 1.23072696 -0.2723291 2.72964520 -0.231469427
## 152 -0.48724019 1.23072696 -0.2723291 2.72964520 -1.253363081
## 153 -0.48724019 1.23072696 3.6647712 2.72964520 -1.811277165
## 154 -0.48724019 1.23072696 -0.2723291 2.72964520 -0.819271765
## 155 -0.48724019 1.23072696 3.6647712 2.72964520 -0.221506675
## 156 -0.48724019 1.23072696 3.6647712 2.72964520 -0.188771920
## 157 -0.48724019 1.23072696 -0.2723291 2.72964520 -1.441232109
## 158 -0.48724019 1.23072696 -0.2723291 0.43412107 0.937018999
## 159 -0.48724019 1.23072696 -0.2723291 0.43412107 -0.311171439
## 160 -0.48724019 1.23072696 -0.2723291 2.72964520 0.320751657
## 161 -0.48724019 1.23072696 3.6647712 0.43412107 -0.049293399
## 162 -0.48724019 1.23072696 -0.2723291 0.43412107 1.714113616
## 163 -0.48724019 1.23072696 3.6647712 0.43412107 2.159590934
## 164 -0.48724019 1.23072696 3.6647712 0.43412107 2.975113306
## 165 -0.48724019 1.23072696 -0.2723291 0.43412107 -0.612900484
## 166 -0.48724019 1.23072696 -0.2723291 0.43412107 -0.261357681
## 167 -0.48724019 1.23072696 -0.2723291 0.43412107 2.340343711
## 168 -0.48724019 1.23072696 -0.2723291 0.43412107 -0.580165729
## 169 -0.48724019 1.23072696 -0.2723291 0.43412107 0.048910866
## 170 -0.48724019 1.23072696 -0.2723291 0.43412107 0.167040633
## 171 -0.48724019 1.23072696 -0.2723291 0.43412107 -0.583012230
## 172 -0.48724019 1.23072696 -0.2723291 0.43412107 -0.575895979
## 173 -0.48724019 -1.03300497 -0.2723291 -0.38570897 -1.014257045
## 174 -0.48724019 -1.03300497 -0.2723291 -0.38570897 0.186966136
## 175 -0.48724019 -1.03300497 -0.2723291 -0.38570897 -0.605784233
## 176 -0.48724019 -1.03300497 -0.2723291 -0.38570897 0.371988664
## 177 -0.48724019 -1.03300497 -0.2723291 -0.38570897 -0.376640949
## 178 -0.48724019 -1.03300497 -0.2723291 -0.38570897 0.043217865
## 179 -0.48724019 -1.03300497 -0.2723291 -0.38570897 0.818889232
## 180 -0.48724019 -1.26477147 -0.2723291 -0.57556435 0.989679257
## 181 -0.48724019 -1.26477147 -0.2723291 -0.57556435 2.106930676
## 182 -0.48724019 -1.26477147 -0.2723291 -0.57556435 -0.200157922
## 183 -0.48724019 -1.26477147 -0.2723291 -0.57556435 1.238748045
## 184 -0.48724019 -1.26477147 -0.2723291 -0.57556435 0.396183918
## 185 -0.48724019 -1.26477147 -0.2723291 -0.57556435 -0.968713038
## 186 -0.48724019 -1.26477147 -0.2723291 -0.57556435 -0.187348670
## 187 -0.48724019 -1.26477147 -0.2723291 -0.57556435 2.200865190
## 188 1.44223095 -1.12192167 -0.2723291 -1.01568364 0.707875715
## 189 1.44223095 -1.12192167 -0.2723291 -1.01568364 0.386221166
## 190 1.44223095 -1.12192167 -0.2723291 -1.01568364 1.281445551
## 191 1.44223095 -1.12192167 -0.2723291 -1.01568364 0.948405001
## 192 1.44223095 -1.12192167 -0.2723291 -1.01568364 0.646675956
## 193 1.44223095 -1.12192167 -0.2723291 -1.01568364 1.271482800
## 194 2.08538800 -1.19626187 -0.2723291 -1.32635608 0.733494219
## 195 2.08538800 -1.19626187 -0.2723291 -1.32635608 0.454537177
## 196 2.94293073 -1.55630166 -0.2723291 -1.14513049 2.263488199
## 197 2.94293073 -1.40179066 -0.2723291 -1.30046671 1.426617073
## 198 2.94293073 -1.40179066 -0.2723291 -1.30046671 1.170432035
## 199 2.94293073 -1.40179066 -0.2723291 -1.30046671 1.408114820
## 200 3.58608778 -1.40907891 -0.2723291 -1.30909650 0.982563006
## 201 3.58608778 -1.40907891 -0.2723291 -1.30909650 1.210283041
## 202 3.05012357 -1.32745046 -0.2723291 -1.20553902 -0.174539418
## 203 3.05012357 -1.32745046 -0.2723291 -1.20553902 1.886326892
## 204 3.58608778 -1.23270315 -0.2723291 -1.19604625 2.232176694
## 205 3.58608778 -1.23270315 -0.2723291 -1.19604625 2.489784983
## 206 -0.48724019 -0.07970124 -0.2723291 -0.56693456 -0.560240226
## 207 -0.48724019 -0.07970124 -0.2723291 -0.56693456 0.058873617
## 208 -0.48724019 -0.07970124 -0.2723291 -0.56693456 -0.713951249
## 209 -0.48724019 -0.07970124 3.6647712 -0.56693456 -0.314017939
## 210 -0.48724019 -0.07970124 3.6647712 -0.56693456 -1.338758093
## 211 -0.48724019 -0.07970124 3.6647712 -0.56693456 -0.462035961
## 212 -0.48724019 -0.07970124 3.6647712 -0.56693456 -1.253363081
## 213 -0.48724019 -0.07970124 3.6647712 -0.56693456 -0.679793244
## 214 -0.48724019 -0.07970124 -0.2723291 -0.56693456 0.128612878
## 215 -0.48724019 -0.07970124 -0.2723291 -0.56693456 -1.241977079
## 216 -0.48724019 -0.07970124 -0.2723291 -0.56693456 -0.146074414
## 217 -0.48724019 0.40132357 3.6647712 -0.04051738 -0.564509977
## 218 -0.48724019 0.40132357 -0.2723291 -0.04051738 0.508620685
## 219 -0.48724019 0.40132357 3.6647712 -0.04051738 -0.474845213
## 220 -0.48724019 0.40132357 3.6647712 -0.04051738 0.125766377
## 221 -0.48724019 -0.71961001 3.6647712 -0.41159834 0.948405001
## 222 -0.48724019 -0.71961001 3.6647712 -0.41159834 -0.171692918
## 223 -0.48724019 -0.71961001 3.6647712 -0.41159834 0.845930986
## 224 -0.48724019 -0.71961001 -0.2723291 -0.41159834 0.474462680
## 225 -0.48724019 -0.71961001 -0.2723291 -0.43748771 2.819979033
## 226 -0.48724019 -0.71961001 -0.2723291 -0.43748771 3.473250881
## 227 -0.48724019 -0.71961001 -0.2723291 -0.43748771 2.498324485
## 228 -0.48724019 -0.71961001 -0.2723291 -0.43748771 1.250134047
## 229 -0.48724019 -0.71961001 -0.2723291 -0.43748771 1.994493909
## 230 -0.48724019 -0.71961001 -0.2723291 -0.43748771 0.380528166
## 231 -0.48724019 -0.71961001 -0.2723291 -0.43748771 -0.432147707
## 232 -0.48724019 -0.71961001 -0.2723291 -0.43748771 1.604523350
## 233 -0.48724019 -0.71961001 -0.2723291 -0.41159834 2.921029798
## 234 -0.48724019 -0.71961001 -0.2723291 -0.41159834 2.792937279
## 235 -0.48724019 -0.71961001 3.6647712 -0.41159834 0.628173703
## 236 -0.48724019 -0.71961001 -0.2723291 -0.41159834 -0.282706434
## 237 -0.48724019 -0.71961001 3.6647712 -0.41159834 0.492964932
## 238 -0.48724019 -0.71961001 -0.2723291 -0.41159834 1.527667838
## 239 0.79907391 -0.90473168 -0.2723291 -1.09335175 0.279477400
## 240 0.79907391 -0.90473168 -0.2723291 -1.09335175 0.457383677
## 241 0.79907391 -0.90473168 -0.2723291 -1.09335175 0.871549489
## 242 0.79907391 -0.90473168 -0.2723291 -1.09335175 -0.269897182
## 243 0.79907391 -0.90473168 -0.2723291 -1.09335175 0.104417624
## 244 0.79907391 -0.90473168 -0.2723291 -1.09335175 0.154231381
## 245 0.45605682 -0.76917014 -0.2723291 -1.06746238 -0.984368790
## 246 0.45605682 -0.76917014 -0.2723291 -1.06746238 -0.967289788
## 247 0.45605682 -0.76917014 -0.2723291 -1.06746238 -0.251394930
## 248 0.45605682 -0.76917014 -0.2723291 -1.06746238 -0.083451404
## 249 0.45605682 -0.76917014 -0.2723291 -1.06746238 0.211161390
## 250 0.45605682 -0.76917014 -0.2723291 -1.06746238 0.616787701
## 251 0.45605682 -0.76917014 -0.2723291 -1.06746238 0.288016902
## 252 0.45605682 -0.76917014 -0.2723291 -1.06746238 0.218277641
## 253 0.45605682 -0.76917014 -0.2723291 -1.06746238 0.956944502
## 254 0.45605682 -0.76917014 -0.2723291 -1.06746238 2.810016281
## 255 2.94293073 -1.09276866 -0.2723291 -1.40402419 -0.251394930
## 256 2.94293073 -1.09276866 -0.2723291 -1.40402419 -0.581588979
## 257 3.37170210 -1.07673449 -0.2723291 -1.38676461 1.664299859
## 258 0.37030254 -1.04466617 -0.2723291 0.79657225 3.443362627
## 259 0.37030254 -1.04466617 -0.2723291 0.79657225 1.492086583
## 260 0.37030254 -1.04466617 -0.2723291 0.79657225 0.793270728
## 261 0.37030254 -1.04466617 -0.2723291 0.79657225 1.307064055
## 262 0.37030254 -1.04466617 -0.2723291 0.79657225 1.758234373
## 263 0.37030254 -1.04466617 -0.2723291 0.79657225 3.007848061
## 264 0.37030254 -1.04466617 -0.2723291 0.79657225 1.483547082
## 265 0.37030254 -1.04466617 -0.2723291 0.79657225 1.311333806
## 266 0.37030254 -1.04466617 -0.2723291 0.79657225 -1.031336047
## 267 0.37030254 -1.04466617 -0.2723291 0.79657225 1.038069765
## 268 0.37030254 -1.04466617 -0.2723291 0.17522737 2.864099790
## 269 0.37030254 -1.04466617 -0.2723291 0.17522737 1.687071862
## 270 0.37030254 -0.60882854 3.6647712 -0.78267931 -0.518965970
## 271 0.37030254 -0.60882854 -0.2723291 -0.78267931 -0.610053984
## 272 0.37030254 -0.60882854 -0.2723291 -0.78267931 -0.063525901
## 273 0.37030254 -0.60882854 -0.2723291 -0.78267931 0.360602663
## 274 0.37030254 -0.60882854 3.6647712 -0.78267931 2.001610160
## 275 1.22784527 -0.68899934 3.6647712 -0.92938574 0.673717710
## 276 1.22784527 -0.68899934 -0.2723291 -0.92938574 0.810349730
## 277 1.22784527 -0.68899934 3.6647712 -0.92938574 1.398152069
## 278 1.22784527 -0.68899934 3.6647712 -0.92938574 0.770498724
## 279 1.22784527 -0.68899934 -0.2723291 -0.92938574 0.280900651
## 280 0.37030254 -1.13795583 -0.2723291 -0.96476788 0.750573221
## 281 0.37030254 -1.13795583 -0.2723291 -0.96476788 2.185209437
## 282 0.37030254 -1.13795583 -0.2723291 -0.96476788 0.972600255
## 283 0.37030254 -1.13795583 3.6647712 -0.96476788 1.936140650
## 284 3.37170210 -1.44697784 3.6647712 -1.32635608 2.331804209
## 285 3.37170210 -1.19043127 -0.2723291 -1.33498587 1.143390280
## 286 1.87100232 -1.29538214 -0.2723291 -1.42991356 0.239626394
## 287 2.94293073 -1.36680703 -0.2723291 -1.46443272 -0.077758404
## 288 1.76380948 -0.84788329 -0.2723291 -1.29183692 -0.107646658
## 289 1.76380948 -0.84788329 -0.2723291 -1.29183692 0.043217865
## 290 1.76380948 -0.84788329 -0.2723291 -1.29183692 0.399030418
## 291 2.94293073 -0.90181638 -0.2723291 -1.24005818 0.820312482
## 292 2.94293073 -0.90181638 -0.2723291 -1.24005818 1.228785293
## 293 2.94293073 -0.90181638 -0.2723291 -1.24005818 0.491541682
## 294 -0.48724019 0.40569652 -0.2723291 -1.01568364 -0.224353176
## 295 -0.48724019 0.40569652 -0.2723291 -1.01568364 -0.392296701
## 296 -0.48724019 0.40569652 -0.2723291 -1.01568364 0.559857693
## 297 -0.48724019 0.40569652 -0.2723291 -1.01568364 0.376258415
## 298 -0.48724019 0.40569652 -0.2723291 -1.01568364 -0.703988498
## 299 2.51415937 -1.29683979 -0.2723291 -1.33498587 0.085915371
## 300 2.51415937 -1.29683979 -0.2723291 -1.33498587 1.076497520
## 301 2.51415937 -1.29683979 -0.2723291 -1.33498587 0.834544984
## 302 0.97058245 -0.73564417 -0.2723291 -1.05020280 0.434611674
## 303 0.97058245 -0.73564417 -0.2723291 -1.05020280 0.299402903
## 304 0.97058245 -0.73564417 -0.2723291 -1.05020280 0.992525758
## 305 0.92770532 -1.30558569 -0.2723291 -0.71364099 1.354031312
## 306 0.92770532 -1.30558569 -0.2723291 -0.71364099 0.471616179
## 307 0.92770532 -1.30558569 -0.2723291 -0.71364099 1.615909352
## 308 0.92770532 -1.30558569 -0.2723291 -0.71364099 0.803233479
## 309 -0.48724019 -0.18027916 -0.2723291 -0.09229612 0.498657933
## 310 -0.48724019 -0.18027916 -0.2723291 -0.09229612 -0.444956959
## 311 -0.48724019 -0.18027916 -0.2723291 -0.09229612 -1.866783923
## 312 -0.48724019 -0.18027916 -0.2723291 -0.09229612 -0.231469427
## 313 -0.48724019 -0.18027916 -0.2723291 -0.09229612 -0.372371198
## 314 -0.48724019 -0.18027916 -0.2723291 -0.09229612 -0.026521396
## 315 -0.48724019 -0.18027916 -0.2723291 -0.09229612 0.401876919
## 316 -0.48724019 -0.18027916 -0.2723291 -0.09229612 -0.824964766
## 317 -0.48724019 -0.18027916 -0.2723291 -0.09229612 -0.527505471
## 318 -0.48724019 -0.18027916 -0.2723291 -0.09229612 -0.715374500
## 319 -0.48724019 -0.18027916 -0.2723291 -0.09229612 0.138575629
## 320 -0.48724019 -0.18027916 -0.2723291 -0.09229612 -0.244278679
## 321 -0.48724019 -0.54760720 -0.2723291 -0.53241540 0.201198639
## 322 -0.48724019 -0.54760720 -0.2723291 -0.53241540 0.130036128
## 323 -0.48724019 -0.54760720 -0.2723291 -0.53241540 -0.346752694
## 324 -0.48724019 -0.54760720 -0.2723291 -0.53241540 -0.820695015
## 325 -0.48724019 -0.54760720 -0.2723291 -0.53241540 0.185542886
## 326 -0.48724019 -0.54760720 -0.2723291 -0.53241540 0.208314890
## 327 -0.48724019 -0.54760720 -0.2723291 -0.53241540 0.038948114
## 328 -0.48724019 -0.54760720 -0.2723291 -0.53241540 -0.286976185
## 329 -0.48724019 -1.15107469 -0.2723291 -0.81719847 -0.592974981
## 330 -0.48724019 -1.15107469 -0.2723291 -0.81719847 0.068836369
## 331 -0.48724019 -1.15107469 -0.2723291 -0.81719847 -0.200157922
## 332 1.01345959 -0.74001712 -0.2723291 -1.00791683 -0.823541516
## 333 1.01345959 -0.74001712 -0.2723291 -1.00791683 -0.360985196
## 334 -0.48724019 -0.86683276 -0.2723291 -0.34256002 0.044641115
## 335 -0.48724019 -0.86683276 -0.2723291 -0.34256002 0.036101614
## 336 -0.48724019 -0.86683276 -0.2723291 -0.34256002 -0.352445695
## 337 -0.48724019 -0.86683276 -0.2723291 -0.34256002 -0.591551731
## 338 -0.48724019 -0.86683276 -0.2723291 -0.34256002 -0.554547225
## 339 -0.48724019 -0.86683276 -0.2723291 -0.34256002 -0.321134190
## 340 -0.48724019 -0.86683276 -0.2723291 -0.34256002 -0.426454706
## 341 -0.48724019 -0.86683276 -0.2723291 -0.34256002 -0.450649960
## 342 1.01345959 -1.40179066 -0.2723291 -0.97253469 1.361147563
## 343 -0.48724019 -1.34785757 -0.2723291 -0.31667065 0.363449163
## 344 1.87100232 -1.07236154 -0.2723291 -0.61008351 0.585476196
## 345 1.87100232 -1.07236154 -0.2723291 -0.61008351 0.838814735
## 346 -0.48724019 -0.98344483 -0.2723291 -0.97253469 -0.385180450
## 347 -0.48724019 -0.98344483 -0.2723291 -0.97253469 -0.550277475
## 348 3.15731641 -1.01842846 -0.2723291 -1.08472196 0.329291158
## 349 2.94293073 -1.33036576 -0.2723291 -1.03294322 0.498657933
## 350 1.22784527 -1.44114723 -0.2723291 -1.08472196 0.931325998
## 351 1.22784527 -1.44114723 -0.2723291 -1.08472196 0.292286652
## 352 2.08538800 -1.37701059 -0.2723291 -1.24005818 0.418955921
## 353 2.08538800 -1.37701059 -0.2723291 -1.24005818 -0.570202978
## 354 3.37170210 -1.32890811 -0.2723291 -1.24868797 0.631020203
## 355 2.94293073 -1.34494227 -0.2723291 -1.22279860 -0.884741275
## 356 2.94293073 -1.34494227 -0.2723291 -1.22279860 -0.496193967
## 357 -0.48724019 1.01499462 3.6647712 1.85803641 -0.103376907
## 358 -0.48724019 1.01499462 3.6647712 1.85803641 0.157077882
## 359 -0.48724019 1.01499462 3.6647712 1.85803641 -0.224353176
## 360 -0.48724019 1.01499462 -0.2723291 1.85803641 -0.245701929
## 361 -0.48724019 1.01499462 -0.2723291 1.85803641 0.161347633
## 362 -0.48724019 1.01499462 -0.2723291 1.85803641 -0.047870149
## 363 -0.48724019 1.01499462 -0.2723291 1.85803641 -1.313139590
## 364 -0.48724019 1.01499462 3.6647712 1.85803641 -0.685486245
## 365 -0.48724019 1.01499462 3.6647712 1.40928733 3.551529643
## 366 -0.48724019 1.01499462 -0.2723291 1.40928733 -3.876413226
## 367 -0.48724019 1.01499462 -0.2723291 1.40928733 -1.881016425
## 368 -0.48724019 1.01499462 -0.2723291 0.65849561 -3.446591661
## 369 -0.48724019 1.01499462 -0.2723291 0.65849561 -1.871053674
## 370 -0.48724019 1.01499462 3.6647712 0.65849561 0.566973944
## 371 -0.48724019 1.01499462 3.6647712 0.65849561 1.040916265
## 372 -0.48724019 1.01499462 -0.2723291 0.65849561 -0.097683907
## 373 -0.48724019 1.01499462 3.6647712 0.97779784 -0.583012230
## 374 -0.48724019 1.01499462 -0.2723291 0.97779784 -1.962141687
## 375 -0.48724019 1.01499462 -0.2723291 0.97779784 -3.055197852
## 376 -0.48724019 1.01499462 -0.2723291 1.00368721 1.463621579
## 377 -0.48724019 1.01499462 -0.2723291 1.00368721 0.518583436
## 378 -0.48724019 1.01499462 -0.2723291 1.00368721 0.724954717
## 379 -0.48724019 1.01499462 -0.2723291 1.00368721 0.135729129
## 380 -0.48724019 1.01499462 -0.2723291 1.00368721 -0.087721155
## 381 -0.48724019 1.01499462 -0.2723291 1.00368721 0.972600255
## 382 -0.48724019 1.01499462 -0.2723291 1.00368721 0.370565414
## 383 -0.48724019 1.01499462 -0.2723291 1.25395112 -1.065494052
## 384 -0.48724019 1.01499462 -0.2723291 1.25395112 -1.088266056
## 385 -0.48724019 1.01499462 -0.2723291 1.25395112 -2.727850303
## 386 -0.48724019 1.01499462 -0.2723291 1.25395112 -1.434115858
## 387 -0.48724019 1.01499462 -0.2723291 1.25395112 -2.323647242
## 388 -0.48724019 1.01499462 -0.2723291 1.25395112 -1.828356167
## 389 -0.48724019 1.01499462 -0.2723291 1.25395112 -1.999146193
## 390 -0.48724019 1.01499462 -0.2723291 1.25395112 -1.273288584
## 391 -0.48724019 1.01499462 -0.2723291 1.25395112 -0.813578764
## 392 -0.48724019 1.01499462 -0.2723291 1.25395112 -0.332520192
## 393 -0.48724019 1.01499462 -0.2723291 1.25395112 -1.777119159
## 394 -0.48724019 1.01499462 -0.2723291 1.19354259 -0.130418661
## 395 -0.48724019 1.01499462 -0.2723291 1.19354259 -0.565933227
## 396 -0.48724019 1.01499462 -0.2723291 1.19354259 0.265244898
## 397 -0.48724019 1.01499462 -0.2723291 1.19354259 0.171310384
## 398 -0.48724019 1.01499462 -0.2723291 1.19354259 -0.765188257
## 399 -0.48724019 1.01499462 -0.2723291 1.19354259 -1.183623820
## 400 -0.48724019 1.01499462 -0.2723291 1.19354259 -0.615746985
## 401 -0.48724019 1.01499462 -0.2723291 1.19354259 -0.423608206
## 402 -0.48724019 1.01499462 -0.2723291 1.19354259 0.083068871
## 403 -0.48724019 1.01499462 -0.2723291 1.19354259 0.169887134
## 404 -0.48724019 1.01499462 -0.2723291 1.19354259 -1.331641842
## 405 -0.48724019 1.01499462 -0.2723291 1.19354259 -1.072610303
## 406 -0.48724019 1.01499462 -0.2723291 1.19354259 -0.856276271
## 407 -0.48724019 1.01499462 -0.2723291 0.90012973 -3.055197852
## 408 -0.48724019 1.01499462 -0.2723291 0.90012973 -0.963020037
## 409 -0.48724019 1.01499462 -0.2723291 0.36508275 -0.950210785
## 410 -0.48724019 1.01499462 -0.2723291 0.36508275 0.807503230
## 411 -0.48724019 1.01499462 -0.2723291 0.36508275 -0.750955755
## 412 -0.48724019 1.01499462 -0.2723291 0.36508275 0.529969438
## 413 -0.48724019 1.01499462 -0.2723291 0.36508275 -2.357805247
## 414 -0.48724019 1.01499462 -0.2723291 0.36508275 -1.607752384
## 415 -0.48724019 1.01499462 -0.2723291 1.19354259 -2.512939520
## 416 -0.48724019 1.01499462 -0.2723291 1.07272553 0.212584640
## 417 -0.48724019 1.01499462 -0.2723291 1.07272553 0.707875715
## 418 -0.48724019 1.01499462 -0.2723291 1.07272553 -1.395688102
## 419 -0.48724019 1.01499462 -0.2723291 1.07272553 -0.466305712
## 420 -0.48724019 1.01499462 -0.2723291 1.40928733 0.767652224
## 421 -0.48724019 1.01499462 -0.2723291 1.40928733 0.179849885
## 422 -0.48724019 1.01499462 -0.2723291 1.40928733 -0.396566452
## 423 -0.48724019 1.01499462 -0.2723291 0.51178918 -0.906090028
## 424 -0.48724019 1.01499462 -0.2723291 0.51178918 -0.258511181
## 425 -0.48724019 1.01499462 -0.2723291 0.25289548 -1.024219796
## 426 -0.48724019 1.01499462 -0.2723291 1.07272553 -0.553123975
## 427 -0.48724019 1.01499462 -0.2723291 0.25289548 -0.637095738
## 428 -0.48724019 1.01499462 -0.2723291 1.07272553 -0.117609410
## 429 -0.48724019 1.01499462 -0.2723291 1.07272553 -0.130418661
## 430 -0.48724019 1.01499462 -0.2723291 1.07272553 0.135729129
## 431 -0.48724019 1.01499462 -0.2723291 0.25289548 0.090185122
## 432 -0.48724019 1.01499462 -0.2723291 0.25289548 0.780461476
## 433 -0.48724019 1.01499462 -0.2723291 0.25289548 0.199775388
## 434 -0.48724019 1.01499462 -0.2723291 1.36613839 0.215431141
## 435 -0.48724019 1.01499462 -0.2723291 1.36613839 -0.109069908
## 436 -0.48724019 1.01499462 -0.2723291 1.59914271 0.490118432
## 437 -0.48724019 1.01499462 -0.2723291 1.59914271 0.251012396
## 438 -0.48724019 1.01499462 -0.2723291 1.59914271 -0.188771920
## 439 -0.48724019 1.01499462 -0.2723291 1.59914271 -0.497617217
## 440 -0.48724019 1.01499462 -0.2723291 1.59914271 -0.935978283
## 441 -0.48724019 1.01499462 -0.2723291 1.59914271 -0.664137492
## 442 -0.48724019 1.01499462 -0.2723291 1.59914271 0.172733634
## 443 -0.48724019 1.01499462 -0.2723291 1.59914271 -0.093414156
## 444 -0.48724019 1.01499462 -0.2723291 1.59914271 0.285170401
## 445 -0.48724019 1.01499462 -0.2723291 1.59914271 -0.612900484
## 446 -0.48724019 1.01499462 -0.2723291 1.59914271 0.248165896
## 447 -0.48724019 1.01499462 -0.2723291 1.59914271 0.080222370
## 448 -0.48724019 1.01499462 -0.2723291 1.59914271 -0.047870149
## 449 -0.48724019 1.01499462 -0.2723291 1.36613839 -0.141804663
## 450 -0.48724019 1.01499462 -0.2723291 1.36613839 0.188389387
## 451 -0.48724019 1.01499462 -0.2723291 1.36613839 0.660908458
## 452 -0.48724019 1.01499462 -0.2723291 1.36613839 0.527122938
## 453 -0.48724019 1.01499462 -0.2723291 1.36613839 0.017599361
## 454 -0.48724019 1.01499462 -0.2723291 1.36613839 1.577481596
## 455 -0.48724019 1.01499462 -0.2723291 1.36613839 0.631020203
## 456 -0.48724019 1.01499462 -0.2723291 1.36613839 0.342100410
## 457 -0.48724019 1.01499462 -0.2723291 1.36613839 -0.439263958
## 458 -0.48724019 1.01499462 -0.2723291 1.36613839 -0.496193967
## 459 -0.48724019 1.01499462 -0.2723291 1.36613839 0.023292362
## 460 -0.48724019 1.01499462 -0.2723291 1.36613839 -0.289822685
## 461 -0.48724019 1.01499462 -0.2723291 1.36613839 0.592592447
## 462 -0.48724019 1.01499462 -0.2723291 1.36613839 0.130036128
## 463 -0.48724019 1.01499462 -0.2723291 1.36613839 0.046064365
## 464 -0.48724019 1.01499462 -0.2723291 1.36613839 0.325021407
## 465 -0.48724019 1.01499462 -0.2723291 0.86561057 -0.107646658
## 466 -0.48724019 1.01499462 -0.2723291 0.86561057 -0.748109254
## 467 -0.48724019 1.01499462 -0.2723291 0.86561057 -0.473421963
## 468 -0.48724019 1.01499462 -0.2723291 0.25289548 -0.400836202
## 469 -0.48724019 1.01499462 -0.2723291 0.21837632 -0.510426469
## 470 -0.48724019 1.01499462 -0.2723291 0.21837632 -0.813578764
## 471 -0.48724019 1.01499462 -0.2723291 0.21837632 -0.167423167
## 472 -0.48724019 1.01499462 -0.2723291 -0.19585359 -0.079181654
## 473 -0.48724019 1.01499462 -0.2723291 0.21837632 0.216854391
## 474 -0.48724019 1.01499462 -0.2723291 0.51178918 0.989679257
## 475 -0.48724019 1.01499462 -0.2723291 0.25289548 -1.220628326
## 476 -0.48724019 1.01499462 -0.2723291 0.25289548 -0.174539418
## 477 -0.48724019 1.01499462 -0.2723291 0.51178918 0.283747151
## 478 -0.48724019 1.01499462 -0.2723291 0.51178918 -1.395688102
## 479 -0.48724019 1.01499462 -0.2723291 0.51178918 -0.141804663
## 480 -0.48724019 1.01499462 -0.2723291 0.51178918 -0.079181654
## 481 -0.48724019 1.01499462 -0.2723291 -0.19585359 -0.060679401
## 482 -0.48724019 1.01499462 -0.2723291 -0.19585359 0.662331708
## 483 -0.48724019 1.01499462 -0.2723291 -0.19585359 1.104962525
## 484 -0.48724019 1.01499462 -0.2723291 -0.19585359 -0.743839504
## 485 -0.48724019 1.01499462 -0.2723291 0.24426569 -0.588705230
## 486 -0.48724019 1.01499462 -0.2723291 0.24426569 0.038948114
## 487 -0.48724019 1.01499462 -0.2723291 0.24426569 -0.242855428
## 488 -0.48724019 1.01499462 -0.2723291 0.24426569 -0.540314723
## 489 -0.48724019 2.42017014 -0.2723291 0.46864023 -1.182200570
## 490 -0.48724019 2.42017014 -0.2723291 0.46864023 -1.239130578
## 491 -0.48724019 2.42017014 -0.2723291 0.46864023 -1.695993897
## 492 -0.48724019 2.42017014 -0.2723291 0.46864023 -0.429301206
## 493 -0.48724019 2.42017014 -0.2723291 0.46864023 -0.429301206
## 494 -0.48724019 -0.21088983 -0.2723291 0.26152527 -0.822118266
## 495 -0.48724019 -0.21088983 -0.2723291 0.26152527 -0.510426469
## 496 -0.48724019 -0.21088983 -0.2723291 0.26152527 -0.874778524
## 497 -0.48724019 -0.21088983 -0.2723291 0.26152527 -1.273288584
## 498 -0.48724019 -0.21088983 -0.2723291 0.26152527 -0.698295497
## 499 -0.48724019 -0.21088983 -0.2723291 0.26152527 -0.378064199
## 500 -0.48724019 -0.21088983 -0.2723291 0.26152527 -1.018526795
## 501 -0.48724019 -0.21088983 -0.2723291 0.26152527 -0.366678197
## 502 -0.48724019 0.11562398 -0.2723291 0.15796779 0.438881424
## 503 -0.48724019 0.11562398 -0.2723291 0.15796779 -0.234315927
## 504 -0.48724019 0.11562398 -0.2723291 0.15796779 0.983986256
## 505 -0.48724019 0.11562398 -0.2723291 0.15796779 0.724954717
## 506 -0.48724019 0.11562398 -0.2723291 0.15796779 -0.362408446
## age dis rad tax ptratio
## 1 -0.119894767 0.1400749840 -0.9818712 -0.66594918 -1.45755797
## 2 0.366803426 0.5566090496 -0.8670245 -0.98635338 -0.30279450
## 3 -0.265548971 0.5566090496 -0.8670245 -0.98635338 -0.30279450
## 4 -0.809087830 1.0766711351 -0.7521778 -1.10502160 0.11292035
## 5 -0.510674339 1.0766711351 -0.7521778 -1.10502160 0.11292035
## 6 -0.350809969 1.0766711351 -0.7521778 -1.10502160 0.11292035
## 7 -0.070159185 0.8384142195 -0.5224844 -0.57694801 -1.50374851
## 8 0.977840575 1.0236248974 -0.5224844 -0.57694801 -1.50374851
## 9 1.116389695 1.0861216287 -0.5224844 -0.57694801 -1.50374851
## 10 0.615481336 1.3283202075 -0.5224844 -0.57694801 -1.50374851
## 11 0.913894827 1.2117799501 -0.5224844 -0.57694801 -1.50374851
## 12 0.508905089 1.1547920492 -0.5224844 -0.57694801 -1.50374851
## 13 -1.050660656 0.7863652700 -0.5224844 -0.57694801 -1.50374851
## 14 -0.240681180 0.4333252240 -0.6373311 -0.60068166 1.17530274
## 15 0.565745754 0.3166899868 -0.6373311 -0.60068166 1.17530274
## 16 -0.428965883 0.3341187865 -0.6373311 -0.60068166 1.17530274
## 17 -1.395257187 0.3341187865 -0.6373311 -0.60068166 1.17530274
## 18 0.466274590 0.2198105553 -0.6373311 -0.60068166 1.17530274
## 19 -1.135921653 0.0006920764 -0.6373311 -0.60068166 1.17530274
## 20 0.032864520 0.0006920764 -0.6373311 -0.60068166 1.17530274
## 21 1.048891406 0.0013569352 -0.6373311 -0.60068166 1.17530274
## 22 0.732715207 0.1031753182 -0.6373311 -0.60068166 1.17530274
## 23 0.821528746 0.0863638874 -0.6373311 -0.60068166 1.17530274
## 24 1.116389695 0.1425444597 -0.6373311 -0.60068166 1.17530274
## 25 0.906789743 0.2871037683 -0.6373311 -0.60068166 1.17530274
## 26 0.608376252 0.3132232229 -0.6373311 -0.60068166 1.17530274
## 27 0.771793164 0.4212152950 -0.6373311 -0.60068166 1.17530274
## 28 0.718505041 0.3126533438 -0.6373311 -0.60068166 1.17530274
## 29 0.917447368 0.3132707128 -0.6373311 -0.60068166 1.17530274
## 30 0.665216917 0.2108349609 -0.6373311 -0.60068166 1.17530274
## 31 0.906789743 0.2079855659 -0.6373311 -0.60068166 1.17530274
## 32 1.116389695 0.1804414138 -0.6373311 -0.60068166 1.17530274
## 33 0.476932215 0.0925850666 -0.6373311 -0.60068166 1.17530274
## 34 0.938762618 -0.0037244859 -0.6373311 -0.60068166 1.17530274
## 35 1.006260907 -0.0167367233 -0.6373311 -0.60068166 1.17530274
## 36 -0.013318520 -0.2064589434 -0.5224844 -0.76681717 0.34387304
## 37 -0.254891346 -0.1981007179 -0.5224844 -0.76681717 0.34387304
## 38 -0.961847117 0.0660856927 -0.5224844 -0.76681717 0.34387304
## 39 -1.363284313 0.0248169544 -0.5224844 -0.76681717 0.34387304
## 40 -1.661697804 0.7627152911 -0.7521778 -0.92701927 -0.07184181
## 41 -1.874850298 0.7627152911 -0.7521778 -0.92701927 -0.07184181
## 42 -2.333128159 0.9145880470 -0.7521778 -1.03975408 -0.25660396
## 43 -2.201684121 0.9145880470 -0.7521778 -1.03975408 -0.25660396
## 44 -2.205236663 0.9145880470 -0.7521778 -1.03975408 -0.25660396
## 45 -1.015135240 0.9145880470 -0.7521778 -1.03975408 -0.25660396
## 46 -1.235392817 0.6199131095 -0.7521778 -1.03975408 -0.25660396
## 47 -1.253155525 0.6199131095 -0.7521778 -1.03975408 -0.25660396
## 48 0.601271169 0.8996287230 -0.7521778 -1.03975408 -0.25660396
## 49 0.949420242 0.9853955139 -0.7521778 -1.03975408 -0.25660396
## 50 -0.233576097 1.0887810641 -0.7521778 -1.03975408 -0.25660396
## 51 -0.812640371 1.4340327636 -0.6373311 -0.98041997 -0.76469989
## 52 -0.198050682 1.4340327636 -0.6373311 -0.98041997 -0.76469989
## 53 -1.686565595 1.4340327636 -0.6373311 -0.98041997 -0.76469989
## 54 -1.675907970 1.4340327636 -0.6373311 -0.98041997 -0.76469989
## 55 -0.745142082 1.6738568465 -0.7521778 0.36053095 1.22149328
## 56 -1.658145262 2.3277455193 -0.5224844 -1.08128796 -0.25660396
## 57 -1.167894527 2.5609210138 -0.8670245 -0.56508119 -0.53374719
## 58 -0.997372532 2.1511780064 -0.5224844 -0.90328562 -1.54993904
## 59 -1.398809729 1.9089794276 -0.1779443 -0.73715011 0.57482574
## 60 -0.759352248 1.4897384367 -0.1779443 -0.73715011 0.57482574
## 61 -0.084369352 1.6290738544 -0.1779443 -0.73715011 0.57482574
## 62 0.881921953 1.4358373805 -0.1779443 -0.73715011 0.57482574
## 63 -0.027528687 1.6291213443 -0.1779443 -0.73715011 0.57482574
## 64 -0.894348827 1.9878601804 -0.1779443 -0.73715011 0.57482574
## 65 -0.322389636 2.5776849546 -0.7521778 -1.14062207 0.06672981
## 66 -1.803799466 1.3375332514 -0.6373311 -0.42267932 -1.08803366
## 67 -1.331311439 1.3375332514 -0.6373311 -0.42267932 -1.08803366
## 68 -1.675907970 1.2836321952 -0.6373311 -0.37521203 0.20530143
## 69 -1.128816570 1.2836321952 -0.6373311 -0.37521203 0.20530143
## 70 -1.263813149 1.2836321952 -0.6373311 -0.37521203 0.20530143
## 71 -2.201684121 0.7086717651 -0.6373311 -0.61254848 0.34387304
## 72 -1.814457091 0.7086717651 -0.6373311 -0.61254848 0.34387304
## 73 -2.159053622 0.7086717651 -0.6373311 -0.61254848 0.34387304
## 74 -2.215894287 0.7086717651 -0.6373311 -0.61254848 0.34387304
## 75 -2.222999370 0.2167712006 -0.5224844 -0.06074124 0.11292035
## 76 -0.837508162 0.3360183832 -0.5224844 -0.06074124 0.11292035
## 77 0.210491598 0.1221237952 -0.5224844 -0.06074124 0.11292035
## 78 -0.809087830 0.1403124336 -0.5224844 -0.06074124 0.11292035
## 79 -0.528437047 0.5789293108 -0.5224844 -0.06074124 0.11292035
## 80 -1.135921653 0.3360183832 -0.5224844 -0.06074124 0.11292035
## 81 -1.246050442 0.7625253315 -0.6373311 -0.75495035 0.25149196
## 82 0.064837394 0.7625253315 -0.6373311 -0.75495035 0.25149196
## 83 -1.292233482 0.7625253315 -0.6373311 -0.75495035 0.25149196
## 84 -0.777114956 0.7625253315 -0.6373311 -0.75495035 0.25149196
## 85 -0.730931915 0.4674704746 -0.7521778 -0.95668632 0.02053927
## 86 -0.443176049 0.3051974268 -0.7521778 -0.95668632 0.02053927
## 87 -0.833955621 0.3002109855 -0.7521778 -0.95668632 0.02053927
## 88 -0.418308258 -0.0225304932 -0.7521778 -0.95668632 0.02053927
## 89 0.629691502 -0.1773001341 -0.8670245 -0.82021787 -0.30279450
## 90 -0.194498140 -0.1807194081 -0.8670245 -0.82021787 -0.30279450
## 91 -0.087921893 -0.3337319220 -0.8670245 -0.82021787 -0.30279450
## 92 0.189176348 -0.3338269018 -0.8670245 -0.82021787 -0.30279450
## 93 -0.531989588 -0.0613297558 -0.6373311 -0.82021787 -0.11803234
## 94 -1.409467353 -0.0613297558 -0.6373311 -0.82021787 -0.11803234
## 95 0.309962761 -0.0855021237 -0.6373311 -0.82021787 -0.11803234
## 96 -0.382782843 -0.1423950448 -0.8670245 -0.78461740 -0.21041342
## 97 0.036417061 -0.1423950448 -0.8670245 -0.78461740 -0.21041342
## 98 0.263779721 -0.1423950448 -0.8670245 -0.78461740 -0.21041342
## 99 -1.125264029 -0.1423950448 -0.8670245 -0.78461740 -0.21041342
## 100 -0.215813389 -0.1423950448 -0.8670245 -0.78461740 -0.21041342
## 101 0.402328842 -0.4830877123 -0.5224844 -0.14380900 1.12911220
## 102 0.096810268 -0.4459031069 -0.5224844 -0.14380900 1.12911220
## 103 0.597718628 -0.5130538501 -0.5224844 -0.14380900 1.12911220
## 104 0.668769459 -0.5130538501 -0.5224844 -0.14380900 1.12911220
## 105 0.761135540 -0.6525317376 -0.5224844 -0.14380900 1.12911220
## 106 0.999155824 -0.8016975682 -0.5224844 -0.14380900 1.12911220
## 107 0.828633829 -0.7522605641 -0.5224844 -0.14380900 1.12911220
## 108 0.590613545 -0.7943366310 -0.5224844 -0.14380900 1.12911220
## 109 1.013365990 -0.6468804374 -0.5224844 -0.14380900 1.12911220
## 110 0.803766038 -0.5935967501 -0.5224844 -0.14380900 1.12911220
## 111 -0.503569256 -0.4830877123 -0.5224844 -0.14380900 1.12911220
## 112 0.462722049 -0.5307200994 -0.4076377 0.14099473 -0.30279450
## 113 0.864159245 -0.6846349217 -0.4076377 0.14099473 -0.30279450
## 114 0.952972784 -0.5922195425 -0.4076377 0.14099473 -0.30279450
## 115 0.555088129 -0.7306526517 -0.4076377 0.14099473 -0.30279450
## 116 0.697189791 -0.6325384823 -0.4076377 0.14099473 -0.30279450
## 117 0.139440767 -0.5057404029 -0.4076377 0.14099473 -0.30279450
## 118 0.498247464 -0.4975246471 -0.4076377 0.14099473 -0.30279450
## 119 0.160756016 -0.6256999342 -0.4076377 0.14099473 -0.30279450
## 120 -0.119894767 -0.4919208369 -0.4076377 0.14099473 -0.30279450
## 121 0.039969603 -0.7300827727 -0.8670245 -1.30675758 0.29768250
## 122 0.551535588 -0.7587191929 -0.8670245 -1.30675758 0.29768250
## 123 0.864159245 -0.8111955516 -0.8670245 -1.30675758 0.29768250
## 124 1.009813449 -0.8788686839 -0.8670245 -1.30675758 0.29768250
## 125 0.967182950 -0.8494724251 -0.8670245 -1.30675758 0.29768250
## 126 0.704294875 -0.8558360740 -0.8670245 -1.30675758 0.29768250
## 127 0.960077867 -0.9677698093 -0.8670245 -1.30675758 0.29768250
## 128 0.974288033 -0.9530004450 -0.6373311 0.17066179 1.26768382
## 129 1.073759197 -0.9415078850 -0.6373311 0.17066179 1.26768382
## 130 0.928104993 -0.8620097633 -0.6373311 0.17066179 1.26768382
## 131 1.077311738 -0.7961887377 -0.6373311 0.17066179 1.26768382
## 132 1.034681240 -0.7237666137 -0.6373311 0.17066179 1.26768382
## 133 1.041786323 -0.6969823003 -0.6373311 0.17066179 1.26768382
## 134 0.952972784 -0.6293091680 -0.6373311 0.17066179 1.26768382
## 135 1.059549031 -0.6881491756 -0.6373311 0.17066179 1.26768382
## 136 1.052443947 -0.7998929513 -0.6373311 0.17066179 1.26768382
## 137 0.885474494 -0.8681834525 -0.6373311 0.17066179 1.26768382
## 138 1.059549031 -0.9237941458 -0.6373311 0.17066179 1.26768382
## 139 1.052443947 -1.0098458762 -0.6373311 0.17066179 1.26768382
## 140 1.041786323 -1.0097983862 -0.6373311 0.17066179 1.26768382
## 141 0.889027036 -1.0367726593 -0.6373311 0.17066179 1.26768382
## 142 1.116389695 -1.1186927669 -0.6373311 0.17066179 1.26768382
## 143 1.116389695 -1.1746358896 -0.5224844 -0.03107419 -1.73470120
## 144 1.116389695 -1.1317999841 -0.5224844 -0.03107419 -1.73470120
## 145 1.038233781 -1.1630958396 -0.5224844 -0.03107419 -1.73470120
## 146 1.116389695 -1.1283332201 -0.5224844 -0.03107419 -1.73470120
## 147 1.116389695 -1.0820305506 -0.5224844 -0.03107419 -1.73470120
## 148 0.963630408 -1.1085299245 -0.5224844 -0.03107419 -1.73470120
## 149 0.896132119 -1.0758568614 -0.5224844 -0.03107419 -1.73470120
## 150 0.935210076 -1.0777089681 -0.5224844 -0.03107419 -1.73470120
## 151 1.020471073 -1.0338757744 -0.5224844 -0.03107419 -1.73470120
## 152 1.116389695 -1.0464131126 -0.5224844 -0.03107419 -1.73470120
## 153 0.690084708 -1.0375799879 -0.5224844 -0.03107419 -1.73470120
## 154 1.063101572 -1.0314062987 -0.5224844 -0.03107419 -1.73470120
## 155 0.974288033 -0.9714740229 -0.5224844 -0.03107419 -1.73470120
## 156 0.498247464 -0.9733261297 -0.5224844 -0.03107419 -1.73470120
## 157 0.903237202 -0.9776477121 -0.5224844 -0.03107419 -1.73470120
## 158 1.024023615 -0.9107344185 -0.5224844 -0.03107419 -1.73470120
## 159 1.116389695 -0.9677223194 -0.5224844 -0.03107419 -1.73470120
## 160 1.116389695 -0.9636381865 -0.5224844 -0.03107419 -1.73470120
## 161 0.853501620 -0.9482039634 -0.5224844 -0.03107419 -1.73470120
## 162 0.789555872 -0.8662838558 -0.5224844 -0.03107419 -1.73470120
## 163 1.052443947 -0.8331358935 -0.5224844 -0.03107419 -1.73470120
## 164 0.899684660 -0.7755306237 -0.5224844 -0.03107419 -1.73470120
## 165 0.825081288 -0.6520568384 -0.5224844 -0.03107419 -1.73470120
## 166 0.867711786 -0.7178778639 -0.5224844 -0.03107419 -1.73470120
## 167 0.981393116 -0.8306664178 -0.5224844 -0.03107419 -1.73470120
## 168 0.377461051 -0.6502047316 -0.5224844 -0.03107419 -1.73470120
## 169 0.977840575 -0.8049743725 -0.5224844 -0.03107419 -1.73470120
## 170 0.945867701 -0.7278032567 -0.5224844 -0.03107419 -1.73470120
## 171 0.924552451 -0.6502047316 -0.5224844 -0.03107419 -1.73470120
## 172 1.020471073 -0.6678709809 -0.5224844 -0.03107419 -1.73470120
## 173 0.707847416 -0.5693768922 -0.5224844 -0.66594918 -0.85708096
## 174 0.551535588 -0.5455369536 -0.5224844 -0.66594918 -0.85708096
## 175 0.004444187 -0.5191325596 -0.5224844 -0.66594918 -0.85708096
## 176 -1.260260608 -0.3147359550 -0.5224844 -0.66594918 -0.85708096
## 177 -0.759352248 -0.1140435641 -0.5224844 -0.66594918 -0.85708096
## 178 0.171413641 -0.2267846280 -0.5224844 -0.66594918 -0.85708096
## 179 0.206939056 -0.4177890758 -0.5224844 -0.66594918 -0.85708096
## 180 -0.361467593 -0.4587728745 -0.7521778 -1.27709053 -0.30279450
## 181 0.523115255 -0.5005640019 -0.7521778 -1.27709053 -0.30279450
## 182 -0.226471014 -0.5685220737 -0.7521778 -1.27709053 -0.30279450
## 183 0.839291454 -0.5197499285 -0.7521778 -1.27709053 -0.30279450
## 184 0.960077867 -0.4502246894 -0.7521778 -1.27709053 -0.30279450
## 185 0.754030456 -0.3833113958 -0.7521778 -1.27709053 -0.30279450
## 186 0.007996729 -0.2447358168 -0.7521778 -1.27709053 -0.30279450
## 187 -0.531989588 -0.2829652003 -0.7521778 -1.27709053 -0.30279450
## 188 -0.976057283 -0.0030596271 -0.5224844 -0.06074124 -1.50374851
## 189 -1.402362270 0.3664594203 -0.5224844 -0.06074124 -1.50374851
## 190 -1.054213197 0.3664594203 -0.5224844 -0.06074124 -1.50374851
## 191 -1.672355429 1.2749890302 -0.5224844 -0.06074124 -1.50374851
## 192 -1.341969064 1.2749890302 -0.5224844 -0.06074124 -1.50374851
## 193 -1.501833434 1.2749890302 -0.5224844 -0.06074124 -1.50374851
## 194 -2.084450250 1.1514202651 -0.9818712 -0.84988492 -1.31898635
## 195 -1.768274051 1.1514202651 -0.9818712 -0.84988492 -1.31898635
## 196 -1.299338565 0.8801578569 -0.6373311 -0.90921904 -1.87327282
## 197 -1.224735192 1.6687754254 -0.8670245 -0.47014661 -2.70470251
## 198 -1.135921653 1.6687754254 -0.8670245 -0.47014661 -2.70470251
## 199 -1.075528447 1.6687754254 -0.8670245 -0.47014661 -2.70470251
## 200 -1.892613005 1.8323307009 -0.7521778 -0.03700760 -0.67231881
## 201 -1.942348587 1.8323307009 -0.7521778 -0.03700760 -0.67231881
## 202 -1.071975905 1.1753551835 -0.8670245 -0.35741180 -1.73470120
## 203 -1.878402839 1.1753551835 -0.8670245 -0.35741180 -1.73470120
## 204 -1.256708066 0.6282713349 -0.6373311 -1.09315478 -1.73470120
## 205 -1.302891107 0.6282713349 -0.6373311 -1.09315478 -1.73470120
## 206 -1.643935096 0.0714045634 -0.6373311 -0.77868399 0.06672981
## 207 -0.571067545 0.2658757752 -0.6373311 -0.77868399 0.06672981
## 208 0.146545850 0.2658757752 -0.6373311 -0.77868399 0.06672981
## 209 -0.336599802 0.2109299408 -0.6373311 -0.77868399 0.06672981
## 210 1.116389695 0.0379716616 -0.6373311 -0.77868399 0.06672981
## 211 0.835738912 0.0389689498 -0.6373311 -0.77868399 0.06672981
## 212 0.711399958 -0.0617571650 -0.6373311 -0.77868399 0.06672981
## 213 -0.524884505 -0.0676459148 -0.6373311 -0.77868399 0.06672981
## 214 -1.288680940 0.0714045634 -0.6373311 -0.77868399 0.06672981
## 215 -2.088002791 -0.0985618510 -0.6373311 -0.77868399 0.06672981
## 216 -0.929874243 0.0714045634 -0.6373311 -0.77868399 0.06672981
## 217 -0.446728591 -0.3243289184 -0.5224844 -0.78461740 -0.94946204
## 218 0.587061003 -0.1775850736 -0.5224844 -0.78461740 -0.94946204
## 219 0.896132119 -0.4301364543 -0.5224844 -0.78461740 -0.94946204
## 220 0.846396537 -0.2050342458 -0.5224844 -0.78461740 -0.94946204
## 221 0.707847416 -0.4432436716 -0.1779443 -0.60068166 -0.48755665
## 222 0.807318580 -0.3547699554 -0.1779443 -0.60068166 -0.48755665
## 223 0.324172928 -0.2483450505 -0.1779443 -0.60068166 -0.48755665
## 224 0.434301716 -0.2483450505 -0.1779443 -0.60068166 -0.48755665
## 225 0.345488177 -0.4277144686 -0.1779443 -0.60068166 -0.48755665
## 226 0.512457630 -0.4277144686 -0.1779443 -0.60068166 -0.48755665
## 227 0.636796585 -0.2751293639 -0.1779443 -0.60068166 -0.48755665
## 228 0.402328842 -0.2751293639 -0.1779443 -0.60068166 -0.48755665
## 229 -1.832219799 -0.1994304356 -0.1779443 -0.60068166 -0.48755665
## 230 -1.675907970 -0.1994304356 -0.1779443 -0.60068166 -0.48755665
## 231 -0.016871062 -0.0586703204 -0.1779443 -0.60068166 -0.48755665
## 232 0.295752595 -0.0586703204 -0.1779443 -0.60068166 -0.48755665
## 233 0.167861099 0.0205903518 -0.1779443 -0.60068166 -0.48755665
## 234 0.064837394 -0.0679783442 -0.1779443 -0.60068166 -0.48755665
## 235 -0.073711727 -0.0679783442 -0.1779443 -0.60068166 -0.48755665
## 236 -0.251338805 -0.0679783442 -0.1779443 -0.60068166 -0.48755665
## 237 0.281542429 0.1676191361 -0.1779443 -0.60068166 -0.48755665
## 238 0.107467893 0.1676191361 -0.1779443 -0.60068166 -0.48755665
## 239 -1.778931675 1.1373157596 -0.4076377 -0.64221553 -0.85708096
## 240 -0.936979326 1.1373157596 -0.4076377 -0.64221553 -0.85708096
## 241 -0.507121797 1.2067460189 -0.4076377 -0.64221553 -0.85708096
## 242 -0.123447309 1.2067460189 -0.4076377 -0.64221553 -0.85708096
## 243 -0.556857379 1.5388905012 -0.4076377 -0.64221553 -0.85708096
## 244 -2.159053622 1.5388905012 -0.4076377 -0.64221553 -0.85708096
## 245 0.281542429 1.9755128019 -0.2927910 -0.46421320 0.29768250
## 246 0.057732311 1.9755128019 -0.2927910 -0.46421320 0.29768250
## 247 -1.196314860 2.0232876588 -0.2927910 -0.46421320 0.29768250
## 248 0.377461051 2.0232876588 -0.2927910 -0.46421320 0.29768250
## 249 -0.691853958 1.9145357480 -0.2927910 -0.46421320 0.29768250
## 250 -1.814457091 1.9145357480 -0.2927910 -0.46421320 0.29768250
## 251 -1.974321461 1.7104240829 -0.2927910 -0.46421320 0.29768250
## 252 -2.119975665 1.7104240829 -0.2927910 -0.46421320 0.29768250
## 253 -2.194579038 2.4275218358 -0.2927910 -0.46421320 0.29768250
## 254 -2.137738373 2.4275218358 -0.2927910 -0.46421320 0.29768250
## 255 -1.299338565 2.5764502168 -0.9818712 -0.55321437 -0.94946204
## 256 -1.757616426 2.5764502168 -0.9818712 -0.55321437 -0.94946204
## 257 -1.221182651 1.2067460189 -0.7521778 -0.97448656 -1.18041474
## 258 0.651006751 -0.9469692255 -0.5224844 -0.85581834 -2.51994036
## 259 1.116389695 -0.9025186628 -0.5224844 -0.85581834 -2.51994036
## 260 1.116389695 -0.8473828687 -0.5224844 -0.85581834 -2.51994036
## 261 0.469827132 -0.7992280924 -0.5224844 -0.85581834 -2.51994036
## 262 0.739820290 -0.7860733853 -0.5224844 -0.85581834 -2.51994036
## 263 0.814423663 -0.7154558781 -0.5224844 -0.85581834 -2.51994036
## 264 0.920999910 -0.8150422349 -0.5224844 -0.85581834 -2.51994036
## 265 0.817976204 -0.8856597421 -0.5224844 -0.85581834 -2.51994036
## 266 -0.205155765 -0.8588754287 -0.5224844 -0.85581834 -2.51994036
## 267 0.569298295 -0.7893501896 -0.5224844 -0.85581834 -2.51994036
## 268 -0.055949019 -0.6522467981 -0.5224844 -0.85581834 -2.51994036
## 269 -0.567515004 -0.4383522101 -0.5224844 -0.85581834 -2.51994036
## 270 -0.251338805 0.0581548764 -0.7521778 -1.09908819 0.06672981
## 271 -0.940531867 0.3010658040 -0.7521778 -1.09908819 0.06672981
## 272 -1.857087590 0.3010658040 -0.7521778 -1.09908819 0.06672981
## 273 -0.350809969 0.0581548764 -0.7521778 -1.09908819 0.06672981
## 274 -0.595935336 0.2713846056 -0.7521778 -1.09908819 0.06672981
## 275 -1.267365691 0.1341862342 -0.6373311 -0.91515245 -0.39517558
## 276 -0.915664077 0.2242746075 -0.6373311 -0.91515245 -0.39517558
## 277 -0.695406500 0.4711746882 -0.6373311 -0.91515245 -0.39517558
## 278 -1.455650394 0.5070770657 -0.6373311 -0.91515245 -0.39517558
## 279 -1.295786023 0.1639624124 -0.6373311 -0.91515245 -0.39517558
## 280 -1.292233482 0.1451564051 -0.5224844 -1.14062207 -1.64232012
## 281 -0.144762558 0.4272465145 -0.5224844 -1.14062207 -1.64232012
## 282 -1.114606404 0.6884410603 -0.5224844 -1.14062207 -1.64232012
## 283 -0.670538709 0.6728643674 -0.5224844 -1.14062207 -1.64232012
## 284 -1.555121557 0.9925190015 -0.9818712 -1.24742347 -2.24279713
## 285 -1.697223220 1.6679680968 -0.9818712 -0.73121670 -1.45755797
## 286 -1.302891107 1.6679680968 -0.9818712 -0.64221553 -1.45755797
## 287 -1.317101273 2.5141909351 -0.9818712 -0.99228679 -0.11803234
## 288 -1.324206356 1.6726695986 -0.4076377 -0.68374941 -0.85708096
## 289 -0.816192913 1.6726695986 -0.4076377 -0.68374941 -0.85708096
## 290 -1.622619847 1.6726695986 -0.4076377 -0.68374941 -0.85708096
## 291 -1.444992769 0.6276539660 -0.6373311 -0.96855315 0.34387304
## 292 -1.452097852 0.6276539660 -0.6373311 -0.96855315 0.34387304
## 293 -1.604857139 0.6276539660 -0.6373311 -0.96855315 0.34387304
## 294 -1.782484217 0.8109650472 -0.6373311 -0.70748306 -1.13422420
## 295 -0.933426784 0.8109650472 -0.6373311 -0.70748306 -1.13422420
## 296 -1.331311439 1.0283263992 -0.6373311 -0.70748306 -1.13422420
## 297 -0.624355669 1.0283263992 -0.6373311 -0.70748306 -1.13422420
## 298 -0.375677759 1.1991001422 -0.6373311 -0.70748306 -1.13422420
## 299 -1.722091011 1.9151531169 -0.5224844 -0.29807769 -1.68851066
## 300 -2.080897708 1.9151531169 -0.5224844 -0.29807769 -1.68851066
## 301 -0.752247165 1.9151531169 -0.5224844 -0.29807769 -1.68851066
## 302 -1.000925074 0.8057411563 -0.2927910 -0.47014661 -1.08803366
## 303 -1.782484217 0.8057411563 -0.2927910 -0.47014661 -1.08803366
## 304 -1.807352008 0.8057411563 -0.2927910 -0.47014661 -1.08803366
## 305 -0.976057283 0.1077818401 -0.2927910 -1.10502160 -0.02565127
## 306 -0.372125218 -0.2018524214 -0.2927910 -1.10502160 -0.02565127
## 307 0.118125517 -0.3304551177 -0.2927910 -1.10502160 -0.02565127
## 308 0.061284852 -0.2908010367 -0.2927910 -1.10502160 -0.02565127
## 309 0.494694923 -0.2267846280 -0.6373311 -0.61848189 -0.02565127
## 310 0.288647512 -0.3288879504 -0.6373311 -0.61848189 -0.02565127
## 311 -1.093291155 -0.6058016588 -0.6373311 -0.61848189 -0.02565127
## 312 -0.560409921 -0.5483863487 -0.6373311 -0.61848189 -0.02565127
## 313 0.775345706 -0.4563983787 -0.6373311 -0.61848189 -0.02565127
## 314 0.505352547 -0.2527616128 -0.6373311 -0.61848189 -0.02565127
## 315 0.665216917 -0.0915333432 -0.6373311 -0.61848189 -0.02565127
## 316 0.324172928 0.0712146037 -0.6373311 -0.61848189 -0.02565127
## 317 0.519562713 0.0966691995 -0.6373311 -0.61848189 -0.02565127
## 318 0.111020434 0.1123883621 -0.6373311 -0.61848189 -0.02565127
## 319 -0.048843936 -0.1246813056 -0.6373311 -0.61848189 -0.02565127
## 320 -0.347257427 0.0982363667 -0.6373311 -0.61848189 -0.02565127
## 321 -0.578172628 0.3539695720 -0.5224844 -0.71934988 0.52863520
## 322 -0.507121797 0.3539695720 -0.5224844 -0.71934988 0.52863520
## 323 -0.663433626 0.4397838527 -0.5224844 -0.71934988 0.52863520
## 324 0.203386515 0.4397838527 -0.5224844 -0.71934988 0.52863520
## 325 -1.011582699 0.4397838527 -0.5224844 -0.71934988 0.52863520
## 326 -1.913928255 0.7697437989 -0.5224844 -0.71934988 0.52863520
## 327 -1.409467353 0.7697437989 -0.5224844 -0.71934988 0.52863520
## 328 -0.883691203 0.7697437989 -0.5224844 -0.71934988 0.52863520
## 329 -1.519596142 0.6741465952 -0.6373311 0.12912791 -0.71850935
## 330 -1.825114716 0.6741465952 -0.6373311 0.12912791 -0.71850935
## 331 -1.292233482 0.9871051509 -0.6373311 0.12912791 -0.71850935
## 332 -1.427230061 1.3514003073 -0.9818712 -0.61848189 -0.71850935
## 333 -1.608409681 1.3514003073 -0.9818712 -0.61848189 -0.71850935
## 334 -1.082633530 1.2648261879 -0.5224844 -1.09315478 0.80577843
## 335 -1.068423364 1.2648261879 -0.5224844 -1.09315478 0.80577843
## 336 -1.210525026 1.0401513886 -0.5224844 -1.09315478 0.80577843
## 337 -0.791325122 0.6819824315 -0.5224844 -1.09315478 0.80577843
## 338 -0.318837095 0.8642962245 -0.5224844 -1.09315478 0.80577843
## 339 -1.111053862 0.4830471675 -0.5224844 -1.09315478 0.80577843
## 340 -0.823297996 0.4830471675 -0.5224844 -1.09315478 0.80577843
## 341 -0.357915052 0.4830471675 -0.5224844 -1.09315478 0.80577843
## 342 -0.684748875 1.5400302593 -0.9818712 -0.73715011 -1.36517689
## 343 -0.315284553 1.1738829960 -0.9818712 0.08166062 -1.18041474
## 344 -0.432518424 0.9199069177 -0.5224844 -0.22687676 -0.39517558
## 345 -1.437887686 1.2681504821 -0.5224844 -0.22687676 -0.39517558
## 346 -0.713169208 2.0033893834 -0.7521778 -0.33367816 0.15911089
## 347 -0.578172628 2.0033893834 -0.7521778 -0.33367816 0.15911089
## 348 -1.452097852 2.2511442825 -0.6373311 -0.33961157 -0.25660396
## 349 -1.381047021 2.1602960705 -0.6373311 -0.76088376 -0.67231881
## 350 -1.210525026 2.3730983904 -0.9818712 -0.43454615 0.57482574
## 351 -0.858823412 2.3730983904 -0.9818712 -0.43454615 0.57482574
## 352 -1.160789444 3.2840499862 -0.6373311 0.01639310 -0.07184181
## 353 -1.778931675 3.2840499862 -0.6373311 0.01639310 -0.07184181
## 354 -1.153684361 3.9566021965 -0.5224844 -1.31269099 -0.67231881
## 355 -1.658145262 3.2248775491 -0.6373311 -0.44047956 1.63720813
## 356 -1.743406260 3.2248775491 -0.6373311 -0.44047956 1.63720813
## 357 1.024023615 -0.7944316108 1.6596029 1.52941294 0.80577843
## 358 0.796660955 -0.6125452271 1.6596029 1.52941294 0.80577843
## 359 0.526667797 -0.5092546567 1.6596029 1.52941294 0.80577843
## 360 0.452064424 -0.6106931203 1.6596029 1.52941294 0.80577843
## 361 0.690084708 -0.6063715378 1.6596029 1.52941294 0.80577843
## 362 0.800213497 -0.7121315839 1.6596029 1.52941294 0.80577843
## 363 0.981393116 -0.8032647354 1.6596029 1.52941294 0.80577843
## 364 0.725610124 -0.8977221812 1.6596029 1.52941294 0.80577843
## 365 0.508905089 -0.8977221812 1.6596029 1.52941294 0.80577843
## 366 0.686532167 -1.0361552904 1.6596029 1.52941294 0.80577843
## 367 0.810871121 -0.9700968153 1.6596029 1.52941294 0.80577843
## 368 1.116389695 -1.0848799457 1.6596029 1.52941294 0.80577843
## 369 1.116389695 -1.1694594886 1.6596029 1.52941294 0.80577843
## 370 1.002708366 -1.1579669285 1.6596029 1.52941294 0.80577843
## 371 1.027576156 -1.2312438711 1.6596029 1.52941294 0.80577843
## 372 1.116389695 -1.2470580136 1.6596029 1.52941294 0.80577843
## 373 0.746925373 -1.2658165310 1.6596029 1.52941294 0.80577843
## 374 1.116389695 -1.2446360278 1.6596029 1.52941294 0.80577843
## 375 1.116389695 -1.2623022771 1.6596029 1.52941294 0.80577843
## 376 1.041786323 -1.1771528552 1.6596029 1.52941294 0.80577843
## 377 0.878369411 -1.1635707388 1.6596029 1.52941294 0.80577843
## 378 1.073759197 -1.1573495596 1.6596029 1.52941294 0.80577843
## 379 0.981393116 -1.1440048928 1.6596029 1.52941294 0.80577843
## 380 1.116389695 -1.1440048928 1.6596029 1.52941294 0.80577843
## 381 0.828633829 -1.1295679579 1.6596029 1.52941294 0.80577843
## 382 1.084416821 -1.0807958128 1.6596029 1.52941294 0.80577843
## 383 1.116389695 -1.0517319833 1.6596029 1.52941294 0.80577843
## 384 1.116389695 -1.0741947142 1.6596029 1.52941294 0.80577843
## 385 0.803766038 -1.1186452769 1.6596029 1.52941294 0.80577843
## 386 1.048891406 -1.1250089259 1.6596029 1.52941294 0.80577843
## 387 1.116389695 -1.1054905698 1.6596029 1.52941294 0.80577843
## 388 0.743372832 -1.0811757321 1.6596029 1.52941294 0.80577843
## 389 1.116389695 -1.0474104008 1.6596029 1.52941294 0.80577843
## 390 1.077311738 -0.9815893753 1.6596029 1.52941294 0.80577843
## 391 1.009813449 -0.8873693792 1.6596029 1.52941294 0.80577843
## 392 0.494694923 -0.7727762084 1.6596029 1.52941294 0.80577843
## 393 1.009813449 -0.9616910999 1.6596029 1.52941294 0.80577843
## 394 0.853501620 -0.9516232374 1.6596029 1.52941294 0.80577843
## 395 0.928104993 -0.9559448199 1.6596029 1.52941294 0.80577843
## 396 1.073759197 -0.9827291333 1.6596029 1.52941294 0.80577843
## 397 0.974288033 -1.0059517029 1.6596029 1.52941294 0.80577843
## 398 1.077311738 -1.0265623271 1.6596029 1.52941294 0.80577843
## 399 1.116389695 -1.0948528283 1.6596029 1.52941294 0.80577843
## 400 0.327725469 -1.0897239172 1.6596029 1.52941294 0.80577843
## 401 1.116389695 -1.0477428302 1.6596029 1.52941294 0.80577843
## 402 1.116389695 -1.0547238481 1.6596029 1.52941294 0.80577843
## 403 1.116389695 -1.0239028917 1.6596029 1.52941294 0.80577843
## 404 0.974288033 -0.9936043244 1.6596029 1.52941294 0.80577843
## 405 0.597718628 -1.0389097056 1.6596029 1.52941294 0.80577843
## 406 1.116389695 -1.1253413553 1.6596029 1.52941294 0.80577843
## 407 1.116389695 -1.2427839210 1.6596029 1.52941294 0.80577843
## 408 1.116389695 -1.1919222195 1.6596029 1.52941294 0.80577843
## 409 1.041786323 -1.1114268095 1.6596029 1.52941294 0.80577843
## 410 1.116389695 -1.1062978984 1.6596029 1.52941294 0.80577843
## 411 1.116389695 -1.1312301050 1.6596029 1.52941294 0.80577843
## 412 1.116389695 -1.0768541496 1.6596029 1.52941294 0.80577843
## 413 1.116389695 -1.0643168114 1.6596029 1.52941294 0.80577843
## 414 1.116389695 -1.0474578907 1.6596029 1.52941294 0.80577843
## 415 1.116389695 -1.0147848276 1.6596029 1.52941294 0.80577843
## 416 1.116389695 -0.9309651233 1.6596029 1.52941294 0.80577843
## 417 0.789555872 -0.9381835908 1.6596029 1.52941294 0.80577843
## 418 0.729162665 -1.0198662487 1.6596029 1.52941294 0.80577843
## 419 1.116389695 -0.9462093868 1.6596029 1.52941294 0.80577843
## 420 0.281542429 -0.9502935197 1.6596029 1.52941294 0.80577843
## 421 1.116389695 -0.9194725633 1.6596029 1.52941294 0.80577843
## 422 0.949420242 -0.9120166463 1.6596029 1.52941294 0.80577843
## 423 0.675874542 -0.8756393696 1.6596029 1.52941294 0.80577843
## 424 0.587061003 -0.8421114879 1.6596029 1.52941294 0.80577843
## 425 0.071942477 -0.8223081923 1.6596029 1.52941294 0.80577843
## 426 0.952972784 -0.8953951752 1.6596029 1.52941294 0.80577843
## 427 -0.315284553 -0.8536040479 1.6596029 1.52941294 0.80577843
## 428 0.359698343 -0.9175729666 1.6596029 1.52941294 0.80577843
## 429 0.338383094 -0.8830477967 1.6596029 1.52941294 0.80577843
## 430 0.960077867 -0.8675660836 1.6596029 1.52941294 0.80577843
## 431 0.622586419 -0.8274371034 1.6596029 1.52941294 0.80577843
## 432 0.913894827 -0.8105781827 1.6596029 1.52941294 0.80577843
## 433 0.221149222 -0.7572944954 1.6596029 1.52941294 0.80577843
## 434 0.686532167 -0.7024911307 1.6596029 1.52941294 0.80577843
## 435 0.938762618 -0.7469416934 1.6596029 1.52941294 0.80577843
## 436 0.924552451 -0.7932443629 1.6596029 1.52941294 0.80577843
## 437 0.878369411 -0.8512295520 1.6596029 1.52941294 0.80577843
## 438 1.116389695 -0.8932106390 1.6596029 1.52941294 0.80577843
## 439 0.686532167 -0.9376612017 1.6596029 1.52941294 0.80577843
## 440 0.899684660 -0.9392758589 1.6596029 1.52941294 0.80577843
## 441 0.846396537 -0.9160057994 1.6596029 1.52941294 0.80577843
## 442 1.016918532 -0.8215483536 1.6596029 1.52941294 0.80577843
## 443 1.116389695 -0.8501847738 1.6596029 1.52941294 0.80577843
## 444 1.116389695 -0.8627221120 1.6596029 1.52941294 0.80577843
## 445 0.995603282 -0.9020437636 1.6596029 1.52941294 0.80577843
## 446 0.931657534 -0.8582105699 1.6596029 1.52941294 0.80577843
## 447 0.988498199 -0.8182715493 1.6596029 1.52941294 0.80577843
## 448 0.995603282 -0.7584342534 1.6596029 1.52941294 0.80577843
## 449 1.070206655 -0.7282306659 1.6596029 1.52941294 0.80577843
## 450 1.055996489 -0.7646079427 1.6596029 1.52941294 0.80577843
## 451 0.853501620 -0.6987869171 1.6596029 1.52941294 0.80577843
## 452 1.052443947 -0.6837801032 1.6596029 1.52941294 0.80577843
## 453 0.825081288 -0.6776064140 1.6596029 1.52941294 0.80577843
## 454 1.091521905 -0.6374774338 1.6596029 1.52941294 0.80577843
## 455 0.906789743 -0.6168668096 1.6596029 1.52941294 0.80577843
## 456 0.636796585 -0.6455032298 1.6596029 1.52941294 0.80577843
## 457 0.686532167 -0.5767378294 1.6596029 1.52941294 0.80577843
## 458 0.416539008 -0.4824228534 1.6596029 1.52941294 0.80577843
## 459 0.537325421 -0.4805707466 1.6596029 1.52941294 0.80577843
## 460 0.562193212 -0.5117241325 1.6596029 1.52941294 0.80577843
## 461 0.761135540 -0.5687120333 1.6596029 1.52941294 0.80577843
## 462 0.704294875 -0.5831489682 1.6596029 1.52941294 0.80577843
## 463 0.512457630 -0.5036983364 1.6596029 1.52941294 0.80577843
## 464 0.757582998 -0.4717851119 1.6596029 1.52941294 0.80577843
## 465 -0.112789684 -0.3949464255 1.6596029 1.52941294 0.80577843
## 466 -0.723826832 -0.3459843207 1.6596029 1.52941294 0.80577843
## 467 0.572850837 -0.4385896596 1.6596029 1.52941294 0.80577843
## 468 0.920999910 -0.5958762661 1.6596029 1.52941294 0.80577843
## 469 0.086152643 -0.4210658801 1.6596029 1.52941294 0.80577843
## 470 -0.421860800 -0.4612898402 1.6596029 1.52941294 0.80577843
## 471 0.547983046 -0.3617034834 1.6596029 1.52941294 0.80577843
## 472 0.786003330 -0.3304076278 1.6596029 1.52941294 0.80577843
## 473 0.228254306 -0.4267171803 1.6596029 1.52941294 0.80577843
## 474 -0.034633770 -0.5993905200 1.6596029 1.52941294 0.80577843
## 475 0.952972784 -0.6483526248 1.6596029 1.52941294 0.80577843
## 476 1.024023615 -0.7546350600 1.6596029 1.52941294 0.80577843
## 477 0.889027036 -0.7074775720 1.6596029 1.52941294 0.80577843
## 478 1.020471073 -0.8046419430 1.6596029 1.52941294 0.80577843
## 479 0.999155824 -0.7714939807 1.6596029 1.52941294 0.80577843
## 480 0.690084708 -0.8756393696 1.6596029 1.52941294 0.80577843
## 481 -0.137657475 -0.1761128861 1.6596029 1.52941294 0.80577843
## 482 0.224701764 -0.2200410597 1.6596029 1.52941294 0.80577843
## 483 0.299305137 -0.1825715149 1.6596029 1.52941294 0.80577843
## 484 -1.004477616 0.1440166471 1.6596029 1.52941294 0.80577843
## 485 -0.947636951 -0.0337381137 1.6596029 1.52941294 0.80577843
## 486 -0.592382795 0.0933923952 1.6596029 1.52941294 0.80577843
## 487 0.398776300 -0.1183176566 1.6596029 1.52941294 0.80577843
## 488 -0.546199754 -0.3052379716 1.6596029 1.52941294 0.80577843
## 489 0.857054162 -0.9375187319 -0.6373311 1.79641644 0.75958789
## 490 1.055996489 -0.9686246278 -0.6373311 1.79641644 0.75958789
## 491 1.045338864 -0.9367114033 -0.6373311 1.79641644 0.75958789
## 492 1.073759197 -0.9151034909 -0.6373311 1.79641644 0.75958789
## 493 0.530220338 -0.8002728706 -0.6373311 1.79641644 0.75958789
## 494 -0.517779422 -0.6711952751 -0.4076377 -0.10227512 0.34387304
## 495 -0.922769160 -0.6711952751 -0.4076377 -0.10227512 0.34387304
## 496 -1.413019895 -0.4732098094 -0.4076377 -0.10227512 0.34387304
## 497 0.153650933 -0.4732098094 -0.4076377 -0.10227512 0.34387304
## 498 0.071942477 -0.4285217972 -0.4076377 -0.10227512 0.34387304
## 499 -0.116342226 -0.6581830377 -0.4076377 -0.10227512 0.34387304
## 500 0.174966182 -0.6625521101 -0.4076377 -0.10227512 0.34387304
## 501 0.395223759 -0.6158695213 -0.4076377 -0.10227512 0.34387304
## 502 0.018654354 -0.6251775451 -0.9818712 -0.80241764 1.17530274
## 503 0.288647512 -0.7159307773 -0.9818712 -0.80241764 1.17530274
## 504 0.796660955 -0.7729186782 -0.9818712 -0.80241764 1.17530274
## 505 0.736267749 -0.6677760011 -0.9818712 -0.80241764 1.17530274
## 506 0.434301716 -0.6126402069 -0.9818712 -0.80241764 1.17530274
## black lstat medv
## 1 0.440615895 -1.074498970 0.159527789
## 2 0.440615895 -0.491952525 -0.101423917
## 3 0.396035074 -1.207532413 1.322937477
## 4 0.415751408 -1.360170785 1.181588636
## 5 0.440615895 -1.025486649 1.486032293
## 6 0.410165113 -1.042290874 0.670558212
## 7 0.426376321 -0.031236706 0.039924923
## 8 0.440615895 0.909799859 0.496590409
## 9 0.328123258 2.419379350 -0.655946292
## 10 0.328999540 0.622727693 -0.394994586
## 11 0.392639483 1.091845624 -0.819041108
## 12 0.440615895 0.086392864 -0.394994586
## 13 0.370513376 0.428078760 -0.090550930
## 14 0.440615895 -0.615183504 -0.231899770
## 15 0.255720500 -0.335113097 -0.471105500
## 16 0.426595391 -0.585776111 -0.286264709
## 17 0.330533033 -0.850442645 0.061670899
## 18 0.329437681 0.282442149 -0.547216415
## 19 -0.741378303 -0.134862757 -0.253645746
## 20 0.375442459 -0.192277190 -0.471105500
## 21 0.217930861 1.171665689 -0.971262937
## 22 0.392749018 0.164812578 -0.318883672
## 23 0.440615895 0.849584722 -0.797295133
## 24 0.414765591 1.012025558 -0.873406047
## 25 0.412465352 0.510699530 -0.753803182
## 26 -0.583319029 0.540106923 -0.938643973
## 27 0.221326452 0.302047077 -0.645073304
## 28 -0.550896614 0.647934029 -0.840787084
## 29 0.342472368 0.020576319 -0.449359525
## 30 0.258020739 -0.094252548 -0.166661844
## 31 0.038293155 1.392921311 -1.069119826
## 32 0.219683424 0.054184768 -0.873406047
## 33 -1.359047220 2.108501199 -1.014754888
## 34 0.022958229 0.797771697 -1.025627875
## 35 -1.186967442 1.076441751 -0.982135924
## 36 0.440615895 -0.416333515 -0.394994586
## 37 0.228774844 -0.174072614 -0.275391721
## 38 0.440615895 -0.543765550 -0.166661844
## 39 0.402607185 -0.353317674 0.235638703
## 40 0.426704926 -1.166922204 0.898890955
## 41 0.426595391 -1.494604580 1.344683452
## 42 0.314759966 -1.094103899 0.442225470
## 43 0.292414788 -0.958269752 0.300876629
## 44 0.413889309 -0.730012370 0.235638703
## 45 0.358354970 -0.434538092 -0.144915868
## 46 0.440615895 -0.342114857 -0.351502635
## 47 0.440615895 0.209623843 -0.275391721
## 48 0.395049257 0.860787538 -0.645073304
## 49 0.440615895 2.542610329 -0.884279035
## 50 0.440615895 0.496696010 -0.340629648
## 51 0.425938180 0.111599201 -0.308010684
## 52 0.408522085 -0.451342316 -0.221026782
## 53 0.440615895 -1.032488409 0.268257666
## 54 0.440615895 -0.591377519 0.094289862
## 55 0.440615895 0.300646725 -0.394994586
## 56 0.429990982 -1.098304955 1.399048391
## 57 0.440615895 -0.963871160 0.235638703
## 58 0.396801820 -1.218735230 0.985874857
## 59 0.372485009 -0.811232788 0.083416874
## 60 0.440615895 -0.480749709 -0.318883672
## 61 0.421009097 0.069588640 -0.416740562
## 62 0.234470674 0.250234052 -0.710311231
## 63 0.440615895 -0.829437365 -0.036185991
## 64 0.426157250 -0.441539852 0.268257666
## 65 0.400526017 -0.644590896 1.138096685
## 66 0.440615895 -1.117909883 0.105162850
## 67 0.440615895 -0.337913801 -0.340629648
## 68 0.433057967 -0.637589136 -0.057931966
## 69 0.440615895 0.061186528 -0.558089402
## 70 0.440615895 -0.540964846 -0.177534831
## 71 0.296358054 -0.830837717 0.181273764
## 72 0.221983663 -0.388326475 -0.090550930
## 73 0.375004318 -0.998879961 0.029051936
## 74 0.224502972 -0.716008850 0.094289862
## 75 0.418927928 -0.822435605 0.170400776
## 76 0.290881295 -0.519959566 -0.123169893
## 77 0.186056121 -0.095652900 -0.275391721
## 78 0.331737920 -0.333712745 -0.188407819
## 79 0.325603949 -0.043839875 -0.144915868
## 80 0.431414939 -0.497553933 -0.242772758
## 81 0.440615895 -1.031088057 0.594447298
## 82 0.426704926 -0.760820115 0.148654801
## 83 0.440615895 -0.830837717 0.246511690
## 84 0.372046868 -0.720209906 0.039924923
## 85 0.440615895 -0.424735627 0.148654801
## 86 0.390229709 -0.857444405 0.442225470
## 87 0.430648193 0.028978431 -0.003567028
## 88 0.421447237 -0.589977167 -0.036185991
## 89 0.440615895 -1.001680665 0.116035838
## 90 0.431414939 -0.973673624 0.670558212
## 91 0.388915287 -0.538164142 0.007305960
## 92 0.403921608 -0.623585616 -0.057931966
## 93 0.419913745 -0.629187024 0.039924923
## 94 0.434372389 -0.902255670 0.268257666
## 95 0.440615895 -0.288901480 -0.210153795
## 96 0.014304949 -0.840640181 0.637939249
## 97 0.385081555 -0.183875078 -0.123169893
## 98 0.440615895 -1.182326077 1.757856987
## 99 0.403702537 -1.271948607 2.312379361
## 100 0.440615895 -0.905056374 1.159842661
## 101 0.417175365 -0.452742668 0.540082359
## 102 0.426157250 -0.697804274 0.431352482
## 103 -3.131326538 -0.283300072 -0.427613550
## 104 0.413998845 0.110198849 -0.351502635
## 105 0.394501581 -0.045240227 -0.264518733
## 106 0.409398367 0.534505515 -0.329756660
## 107 0.427143067 0.841182610 -0.329756660
## 108 0.339733988 0.201221731 -0.231899770
## 109 0.422433054 -0.053642339 -0.297137697
## 110 0.378509444 0.405673128 -0.340629648
## 111 0.403264396 0.048583360 -0.090550930
## 112 0.426266786 -0.349116618 0.029051936
## 113 0.419256534 0.498096362 -0.405867574
## 114 0.440615895 0.621327341 -0.416740562
## 115 0.351235183 -0.308506409 -0.438486537
## 116 -0.128857540 0.435080520 -0.460232513
## 117 0.401183228 -0.085850436 -0.144915868
## 118 0.414436985 -0.329511689 -0.362375623
## 119 -0.197645637 0.380466791 -0.231899770
## 120 0.381466894 0.134004834 -0.351502635
## 121 0.355726125 0.240431588 -0.057931966
## 122 0.229979731 0.226428068 -0.242772758
## 123 0.234580209 0.738956911 -0.221026782
## 124 0.149361834 1.786420232 -0.568962390
## 125 0.248710248 0.689944590 -0.405867574
## 126 0.310488093 0.302047077 -0.123169893
## 127 0.028654058 2.045485358 -0.742930194
## 128 0.388148541 0.635330861 -0.688565255
## 129 0.440615895 0.383267495 -0.492851476
## 130 0.440615895 0.796371345 -0.895152022
## 131 0.420242350 -0.007430722 -0.362375623
## 132 0.440615895 -0.055042691 -0.318883672
## 133 0.318593697 -0.214682823 0.050797911
## 134 0.350687507 0.332854822 -0.449359525
## 135 -1.028689097 0.652135085 -0.753803182
## 136 0.416189548 0.603122764 -0.481978488
## 137 0.236332772 0.594720652 -0.558089402
## 138 0.409726972 0.271239333 -0.590708366
## 139 0.387381794 1.213676250 -1.003881900
## 140 0.440615895 0.813175569 -0.514597451
## 141 0.344005860 1.611376228 -0.927770986
## 142 0.440615895 3.046737061 -0.884279035
## 143 0.440615895 1.983869868 -0.993008912
## 144 0.440615895 1.927855787 -0.753803182
## 145 0.440615895 2.329756820 -1.166976716
## 146 -2.012862748 2.121104367 -0.949516961
## 147 -2.052733556 0.559711851 -0.753803182
## 148 0.383767133 2.363365269 -0.862533059
## 149 0.003460966 2.193922673 -0.514597451
## 150 -0.052840120 1.231880827 -0.775549157
## 151 0.176636095 0.202622083 -0.112296905
## 152 -0.165113687 0.087793216 -0.318883672
## 153 -0.146711775 -0.074647619 -0.786422145
## 154 -1.037561447 0.439281577 -0.340629648
## 155 -0.390537100 0.345457990 -0.601581353
## 156 -2.942816482 0.331454470 -0.753803182
## 157 -2.936025300 0.488293898 -1.025627875
## 158 0.074001626 -1.129112700 2.040554668
## 159 -0.030494942 -0.871447926 0.192146752
## 160 0.083640722 -0.737014131 0.083416874
## 161 -0.194469117 -1.001680665 0.485717421
## 162 0.194490331 -1.529613381 2.986504601
## 163 0.360764744 -1.503006692 2.986504601
## 164 0.348058662 -1.306957408 2.986504601
## 165 0.421009097 -0.141864517 0.018178948
## 166 -1.276238619 -0.398128939 0.268257666
## 167 0.138298780 -1.253744030 2.986504601
## 168 -1.413705278 -0.071846915 0.137781813
## 169 -0.652654802 -0.217483527 0.137781813
## 170 -0.291736362 -0.186675782 -0.025313003
## 171 -0.705231691 0.248833700 -0.558089402
## 172 -0.093587210 -0.087250788 -0.373248611
## 173 0.440615895 0.285242853 0.061670899
## 174 0.425280969 -0.505956045 0.116035838
## 175 0.400416482 -0.421934923 0.007305960
## 176 0.375551994 -1.025486649 0.746669127
## 177 0.400416482 -0.356118378 0.072543887
## 178 0.426376321 -0.891052854 0.224765715
## 179 0.378947585 -0.802830676 0.801034065
## 180 0.440615895 -1.066096858 1.594762170
## 181 0.425938180 -0.713208146 1.877459852
## 182 0.440615895 -0.448541612 1.486032293
## 183 0.410165113 -1.096904603 1.670873085
## 184 0.440615895 -0.976474328 1.083731747
## 185 0.375990135 0.185817859 0.420479494
## 186 0.333380947 0.069588640 0.768415102
## 187 0.393844370 -1.148717628 2.986504601
## 188 0.407426733 -0.836439125 1.029366808
## 189 0.286609423 -1.133313756 0.790161078
## 190 0.440615895 -1.017084537 1.344683452
## 191 0.230089266 -1.057694746 1.573016195
## 192 0.361860096 -1.115109179 0.866271992
## 193 0.370403840 -1.369973249 1.507778268
## 194 0.401949974 -1.067497210 0.931509918
## 195 0.219354818 -1.158520092 0.714050163
## 196 0.411370000 -1.355969729 2.986504601
## 197 0.440615895 -1.200530653 1.170715648
## 198 -0.025894464 -0.566171183 0.844526016
## 199 0.389134357 -0.844841237 1.312064489
## 200 0.440615895 -1.133313756 1.344683452
## 201 0.302601560 -1.148717628 1.127223698
## 202 0.406331382 -0.731412722 0.170400776
## 203 0.423966547 -1.336364800 2.149284545
## 204 0.395487398 -1.238340158 2.823409785
## 205 0.371061052 -1.368572897 2.986504601
## 206 0.440615895 -0.249691623 0.007305960
## 207 0.418380252 -0.235688103 0.203019739
## 208 0.358793111 0.757161488 -0.003567028
## 209 0.269960074 0.281041797 0.203019739
## 210 0.440615895 1.461538560 -0.275391721
## 211 0.400635552 0.646533677 -0.090550930
## 212 0.422433054 1.586169891 -0.351502635
## 213 0.375332924 0.472890025 -0.014440015
## 214 0.319141373 -0.458344076 0.605320286
## 215 -0.084824395 2.366165973 0.126908825
## 216 0.404797889 -0.445740908 0.268257666
## 217 0.395706469 0.120001313 0.083416874
## 218 0.395487398 -0.414933163 0.670558212
## 219 0.440615895 0.737556559 -0.112296905
## 220 0.406002776 -0.301504649 0.050797911
## 221 0.383657598 -0.412132459 0.453098458
## 222 0.422433054 1.233281179 -0.090550930
## 223 0.369308489 -0.381324714 0.540082359
## 224 0.440615895 -0.707606738 0.822780041
## 225 0.310816699 -1.192128541 2.421109239
## 226 0.277408467 -1.123511291 2.986504601
## 227 0.336338397 -1.333564096 1.638254121
## 228 0.168749562 -0.881250390 0.985874857
## 229 0.228227168 -1.222936286 2.627696006
## 230 0.259225626 -1.245341918 0.975001869
## 231 0.237428124 -0.140464165 0.192146752
## 232 0.213220848 -1.036689465 0.996747845
## 233 0.320236725 -1.425987330 2.084046619
## 234 0.244000235 -1.218735230 2.801663810
## 235 0.038621760 -0.644590896 0.703177176
## 236 0.219902494 -0.248291271 0.159527789
## 237 0.348058662 -0.435938444 0.279130654
## 238 0.365803363 -1.109507771 0.975001869
## 239 0.249038854 -0.881250390 0.126908825
## 240 0.296905730 -0.739814835 0.083416874
## 241 0.378728515 -0.178273670 -0.057931966
## 242 0.415641872 -0.035437762 -0.264518733
## 243 0.176088420 -0.200679302 -0.036185991
## 244 0.197557316 -1.045091578 0.126908825
## 245 0.173240505 -0.021434242 -0.536343427
## 246 0.355507055 0.813175569 -0.438486537
## 247 0.367008250 -0.489151821 0.192146752
## 248 0.213220848 -0.350516970 -0.221026782
## 249 0.197557316 -0.438739148 0.213892727
## 250 0.406002776 -0.853243349 0.398733519
## 251 0.433824713 -0.945666583 0.203019739
## 252 0.223407620 -1.269147903 0.246511690
## 253 0.322208358 -1.277550015 0.768415102
## 254 0.440615895 -1.276149663 2.203649484
## 255 0.396692285 -0.851842997 -0.068804954
## 256 0.421775843 -0.476548653 -0.177534831
## 257 0.324946738 -1.336364800 2.334125337
## 258 0.361750561 -1.054894042 2.986504601
## 259 0.291538506 -0.681000049 1.464286318
## 260 0.386176907 -0.805631380 0.822780041
## 261 0.395706469 -0.428936683 1.225080587
## 262 0.347182381 -0.755218707 2.236268447
## 263 0.330642568 -0.944266231 2.856028748
## 264 0.402497650 -0.196478246 0.920636930
## 265 0.341924692 -0.637589136 1.518651256
## 266 0.391325061 -0.308506409 0.029051936
## 267 0.300082251 0.299246373 0.888017967
## 268 0.305230404 -0.730012370 2.986504601
## 269 0.368322672 -1.329363040 2.279760398
## 270 0.379714331 0.139606242 -0.199280807
## 271 0.350249366 0.048583360 -0.155788856
## 272 0.440615895 -0.849042293 0.290003641
## 273 0.419366069 -0.689402161 0.203019739
## 274 0.373470826 -0.850442645 1.377302416
## 275 0.440615895 -1.277550015 1.072858759
## 276 0.440615895 -1.354569377 1.029366808
## 277 0.356821477 -0.924661303 1.159842661
## 278 0.402826256 -1.189327837 1.148969673
## 279 0.440615895 -0.765021171 0.714050163
## 280 0.440615895 -1.092703547 1.366429428
## 281 0.335571651 -1.245341918 2.486347165
## 282 0.389462963 -1.129112700 1.399048391
## 283 0.223407620 -1.350368321 2.551585092
## 284 0.425500039 -1.329363040 2.986504601
## 285 0.416737224 -0.672597937 1.051112783
## 286 0.416737224 -0.619384560 -0.057931966
## 287 -0.165113687 0.038780895 -0.264518733
## 288 0.440615895 -0.772022931 0.072543887
## 289 0.440615895 -0.707606738 -0.025313003
## 290 0.164806295 -0.440139500 0.246511690
## 291 0.440615895 -1.305557056 0.648812237
## 292 0.440615895 -1.273348959 1.605635158
## 293 0.440615895 -1.113708827 0.583574310
## 294 0.440615895 -0.570372239 0.148654801
## 295 0.440615895 -0.315508169 -0.090550930
## 296 0.440615895 -0.893853558 0.659685225
## 297 0.396254144 -0.737014131 0.496590409
## 298 0.440615895 0.446283337 -0.242772758
## 299 0.126688050 -1.075899322 -0.003567028
## 300 0.163272803 -1.108107419 0.703177176
## 301 0.374456642 -0.921860599 0.246511690
## 302 0.428019349 -0.441539852 -0.057931966
## 303 0.295043632 -0.557769070 0.420479494
## 304 0.369746629 -1.091303195 1.148969673
## 305 0.405345565 -0.801430324 1.475159305
## 306 0.401840439 -0.521359918 0.637939249
## 307 0.440615895 -0.865846518 1.181588636
## 308 0.440615895 -0.717409202 0.616193274
## 309 0.440615895 -1.136114460 0.029051936
## 310 0.433386573 -0.375723306 -0.242772758
## 311 -0.068175046 -0.001829314 -0.699438243
## 312 0.440615895 -0.934463767 -0.047058979
## 313 0.434043784 -0.130661701 -0.340629648
## 314 0.402169045 -0.665596177 -0.101423917
## 315 0.427362137 -0.472347596 0.137781813
## 316 0.435358206 -0.161469445 -0.688565255
## 317 0.372704079 0.794970993 -0.514597451
## 318 0.440615895 0.460286857 -0.297137697
## 319 0.422104448 -0.321109577 0.061670899
## 320 0.433277038 0.010773855 -0.166661844
## 321 0.440615895 -0.763620819 0.137781813
## 322 0.440615895 -0.809832436 0.061670899
## 323 0.440615895 -0.693603218 -0.231899770
## 324 0.377414092 -0.127860997 -0.438486537
## 325 0.440615895 -0.914858839 0.268257666
## 326 0.405345565 -1.060495450 0.224765715
## 327 0.440615895 -0.910657783 0.050797911
## 328 0.440615895 0.019175967 -0.036185991
## 329 0.282228015 -0.375723306 -0.351502635
## 330 0.203034075 -0.744015891 0.007305960
## 331 0.130302712 -0.498954285 -0.297137697
## 332 0.409069761 -0.031236706 -0.590708366
## 333 0.061076474 -0.675398641 -0.340629648
## 334 0.361860096 -0.976474328 -0.036185991
## 335 0.358464505 -0.826636661 -0.199280807
## 336 0.440615895 -0.650192305 -0.155788856
## 337 0.440615895 -0.399529291 -0.329756660
## 338 0.417723041 -0.293102536 -0.438486537
## 339 0.432291221 -0.580174703 -0.210153795
## 340 0.440615895 -0.407931403 -0.384121599
## 341 0.440615895 -0.470947244 -0.416740562
## 342 0.416956295 -1.003081017 1.105477722
## 343 0.364598476 -0.560569774 -0.655946292
## 344 0.440615895 -0.766421523 0.148654801
## 345 0.342800973 -1.126311996 0.942382906
## 346 0.317279275 -0.297303592 -0.547216415
## 347 0.086926778 0.002371742 -0.579835378
## 348 0.391653667 -0.881250390 0.061670899
## 349 0.375332924 -0.933063415 0.213892727
## 350 0.363393588 -0.947066935 0.442225470
## 351 0.440615895 -0.934463767 0.039924923
## 352 0.154509988 -1.003081017 0.170400776
## 353 0.390558315 -0.681000049 -0.427613550
## 354 0.304354123 -1.141715868 0.822780041
## 355 0.286171282 -0.644590896 -0.471105500
## 356 0.212125496 -0.991878200 -0.210153795
## 357 0.230636942 0.692745294 -0.514597451
## 358 0.379714331 0.086392864 -0.090550930
## 359 0.424514223 -0.164270149 0.018178948
## 360 0.373142220 0.002371742 0.007305960
## 361 0.195914288 -0.681000049 0.268257666
## 362 -0.065984343 0.215225251 -0.286264709
## 363 0.264154709 -0.344915562 -0.188407819
## 364 -0.039805433 0.278241093 -0.623327329
## 365 -0.023265620 -1.031088057 -0.068804954
## 366 -0.021622592 -0.774823635 0.540082359
## 367 -0.445195159 0.188618563 -0.068804954
## 368 -2.467324237 0.094794977 0.061670899
## 369 0.206429666 -1.315359520 2.986504601
## 370 0.204348498 -1.249542974 2.986504601
## 371 0.387491330 -1.357370081 2.986504601
## 372 0.103795196 -0.437338796 2.986504601
## 373 -0.096325589 -0.528361678 2.986504601
## 374 0.440615895 3.097149734 -0.949516961
## 375 0.440615895 3.545262384 -0.949516961
## 376 0.440615895 0.110198849 -0.819041108
## 377 0.069510683 1.482543841 -0.938643973
## 378 0.440615895 1.202473434 -1.003881900
## 379 0.440615895 1.545559682 -1.025627875
## 380 0.406002776 1.278092444 -1.340944520
## 381 0.440615895 0.638131565 -1.319198544
## 382 0.440615895 1.180067802 -1.264833606
## 383 0.440615895 1.532956514 -1.221341655
## 384 0.440615895 1.667390309 -1.112611777
## 385 -0.775991422 2.517403992 -1.493166348
## 386 0.440615895 2.542610329 -1.667134152
## 387 0.440615895 2.188321265 -1.308325557
## 388 0.440615895 2.707851869 -1.645388177
## 389 0.177950518 2.516003640 -1.340944520
## 390 0.440615895 1.147859705 -1.199595679
## 391 0.413560704 0.624128045 -0.808168120
## 392 0.237756730 0.855186130 0.072543887
## 393 0.440615895 1.824229736 -1.395309459
## 394 0.440615895 0.352459751 -0.949516961
## 395 0.440615895 0.517701290 -1.069119826
## 396 0.386724583 0.625528397 -1.025627875
## 397 0.440615895 0.940607604 -1.090865802
## 398 0.398992524 1.017626966 -1.525785311
## 399 0.440615895 2.511802584 -1.906339882
## 400 -0.202793791 2.424980758 -1.764991042
## 401 0.440615895 1.976868108 -1.841101956
## 402 0.440615895 1.073641047 -1.667134152
## 403 0.212892242 1.072240695 -1.134357753
## 404 0.440615895 0.996621685 -1.547531287
## 405 -0.298089403 2.062289582 -1.525785311
## 406 0.309940417 1.446134688 -1.906339882
## 407 0.148376017 1.496547361 -1.156103728
## 408 -0.269281649 -0.073247267 0.583574310
## 409 -0.460420549 1.925055083 -0.579835378
## 410 -1.942212553 0.998022037 0.540082359
## 411 -3.878356510 -0.356118378 -0.819041108
## 412 -3.522914830 1.199672730 -0.579835378
## 413 -3.591483857 3.041135653 -0.503724464
## 414 -1.595971828 1.040032598 -0.677692268
## 415 -2.939968567 3.406627533 -1.688880128
## 416 -3.608352275 2.296148371 -1.667134152
## 417 -3.670568261 1.839633609 -1.634515189
## 418 -2.511795523 1.958663532 -1.319198544
## 419 -3.726650277 1.115651608 -1.493166348
## 420 -3.376137680 1.412526239 -1.536658299
## 421 -0.415401588 0.331454470 -0.634200317
## 422 -0.401928760 0.426678408 -0.906025010
## 423 -0.713337295 0.202622083 -0.188407819
## 424 -3.879232792 1.489545601 -0.993008912
## 425 -3.866855315 0.631129805 -1.177849704
## 426 -3.822712635 1.643584324 -1.547531287
## 427 -3.636831424 0.425278056 -1.340944520
## 428 -3.700690438 0.261436868 -1.264833606
## 429 -2.847301799 1.241683291 -1.253960618
## 430 -3.241738006 1.600173411 -1.417055434
## 431 -2.992764527 0.698346703 -0.873406047
## 432 -3.015985986 0.985418869 -0.916897998
## 433 -2.833938506 -0.087250788 -0.699438243
## 434 -2.809402625 0.499496714 -0.895152022
## 435 -2.804583076 0.352459751 -1.177849704
## 436 -2.703591634 1.486744897 -0.993008912
## 437 -3.605723431 0.755761136 -1.406182446
## 438 -3.804748865 1.932056843 -1.504039336
## 439 -3.151590547 2.992123331 -1.536658299
## 440 0.440615895 1.432131167 -1.058246839
## 441 0.380919218 1.324304061 -1.308325557
## 442 0.320784401 0.961612885 -0.590708366
## 443 0.427362137 0.551309739 -0.449359525
## 444 0.329218610 0.867789298 -0.775549157
## 445 -1.272295352 1.559563202 -1.275706593
## 446 -3.435177145 1.586169891 -1.166976716
## 447 -0.423507192 0.719351983 -0.829914096
## 448 0.348825409 0.530304459 -1.079992814
## 449 0.440615895 0.766963952 -0.916897998
## 450 -0.574665749 0.932205492 -1.036500863
## 451 -3.903330533 0.670339662 -0.993008912
## 452 -0.015160016 0.710949871 -0.797295133
## 453 0.311254840 0.646533677 -0.699438243
## 454 0.210263398 0.572315020 -0.514597451
## 455 -3.833666154 0.848184370 -0.829914096
## 456 -3.349082489 0.766963952 -0.916897998
## 457 -3.792042783 0.890194931 -1.069119826
## 458 -3.868498343 0.600322060 -0.982135924
## 459 -0.925178346 0.500897066 -0.829914096
## 460 0.440615895 0.286643205 -0.275391721
## 461 -1.111169093 0.527503755 -0.666819280
## 462 0.380700148 0.279641445 -0.525470439
## 463 0.440615895 0.187218211 -0.329756660
## 464 0.406879058 -0.330912041 -0.253645746
## 465 0.440615895 0.079391104 -0.123169893
## 466 -0.243979021 0.206823139 -0.286264709
## 467 -3.665748713 0.629729453 -0.384121599
## 468 -0.278044464 1.213676250 -0.373248611
## 469 0.132164810 0.766963952 -0.373248611
## 470 0.440615895 0.295045317 -0.264518733
## 471 0.440615895 0.509299178 -0.286264709
## 472 0.423418871 0.030378783 -0.318883672
## 473 0.401949974 0.239031236 0.072543887
## 474 0.197228711 -0.139063813 0.790161078
## 475 -0.044844052 0.768364304 -0.949516961
## 476 -0.590548351 1.602974115 -1.003881900
## 477 0.433057967 0.843983314 -0.634200317
## 478 -0.078799960 1.716402630 -1.145230740
## 479 0.252215374 0.752960432 -0.862533059
## 480 0.291867112 0.063987232 -0.123169893
## 481 0.440615895 -0.267896200 0.050797911
## 482 0.398663919 -0.688001809 0.126908825
## 483 0.422871195 -0.790227508 0.268257666
## 484 0.397020891 -0.312707465 -0.079677942
## 485 0.153962312 0.096195329 -0.210153795
## 486 0.349920761 -0.290301832 -0.144915868
## 487 0.394392046 0.325853062 -0.373248611
## 488 0.345539353 -0.168471205 -0.210153795
## 489 0.420790026 0.757161488 -0.797295133
## 490 -0.138277566 1.584769539 -1.688880128
## 491 -0.418906714 2.384370549 -1.569277262
## 492 0.366241503 0.758561840 -0.971262937
## 493 0.440615895 0.097595681 -0.264518733
## 494 0.440615895 -0.090051492 -0.079677942
## 495 0.440615895 0.131204129 0.213892727
## 496 0.401073693 0.692745294 0.061670899
## 497 0.440615895 1.188469914 -0.308010684
## 498 0.440615895 0.202622083 -0.460232513
## 499 0.440615895 0.037380543 -0.144915868
## 500 0.428238419 0.342657286 -0.547216415
## 501 0.440615895 0.234830180 -0.623327329
## 502 0.386834118 -0.417733867 -0.014440015
## 503 0.440615895 -0.500354637 -0.210153795
## 504 0.440615895 -0.982075736 0.148654801
## 505 0.402826256 -0.864446166 -0.057931966
## 506 0.440615895 -0.668396881 -1.156103728
boston_scaled<-data.frame(boston_scaled, crime)
n<-nrow(boston_scaled)
ind<-sample(n, size = n*0.8)
train<-boston_scaled[ind, ]
test<-boston_scaled[-ind, ]
correct_classes<-test$crime
lda.fit<-lda(crime~., data=train)
lda.fit
## Call:
## lda(crime ~ ., data = train)
##
## Prior probabilities of groups:
## low med_low med_high high
## 0.2722772 0.2400990 0.2524752 0.2351485
##
## Group means:
## zn indus chas nox rm
## low 1.0333390 -0.93297036 -0.12916179 -0.8924424 0.44137858
## med_low -0.1360455 -0.25174914 -0.06938576 -0.5372285 -0.14928773
## med_high -0.3686975 0.09798927 0.15226017 0.4070472 0.04159927
## high -0.4872402 1.01726549 -0.06511327 1.0388422 -0.34721712
## age dis rad tax ptratio
## low -0.9028749 0.9339009 -0.6989307 -0.7377435 -0.46110207
## med_low -0.3289820 0.2567362 -0.5556360 -0.4445167 -0.02422269
## med_high 0.3710177 -0.3671523 -0.4267788 -0.3396116 -0.32815401
## high 0.7875739 -0.8325600 1.6366336 1.5129868 0.77903654
## black lstat medv
## low 0.37944248 -0.77353786 0.529604753
## med_low 0.30982524 -0.11648494 -0.009507938
## med_high 0.06910905 0.02974725 0.164964282
## high -0.76317926 0.83687837 -0.667048185
##
## Coefficients of linear discriminants:
## LD1 LD2 LD3
## zn 0.09243423 0.66571366 -0.76537402
## indus 0.01692218 -0.08931198 0.33328572
## chas -0.03780800 -0.05265003 0.07584029
## nox 0.40204335 -0.70391185 -1.53724859
## rm 0.04361795 0.05629763 -0.04717913
## age 0.25790284 -0.26742966 -0.21306498
## dis -0.01270129 -0.12003686 -0.16505629
## rad 3.38286831 0.96987665 -0.24681113
## tax 0.03544286 -0.09253273 0.70753667
## ptratio 0.13160706 0.03888716 -0.22843132
## black -0.11244411 0.05102551 0.15710993
## lstat 0.21020556 -0.21320942 0.43265709
## medv 0.09379916 -0.39166010 -0.22347964
##
## Proportion of trace:
## LD1 LD2 LD3
## 0.9515 0.0355 0.0130
lda.arrows<-function(x, myscale=1, arrow_heads= 0.1, color="red", tex=0.75, choices =c (1,2)){
heads<-coef(x)
arrows(x0 = 0, y0=0, x1 = myscale*heads[,choices[1]],y1=myscale*heads[,choices[2]], col=color, length= arrow_heads)
text(myscale * heads[,choices], labels = row.names(heads),
cex = tex, col=color, pos=3)
}
classes<-as.numeric(train$crime)
plot(lda.fit, dimen=2, col=classes, pch=classes)
lda.arrows(lda.fit, myscale=1)
lda.pred<-predict(lda.fit, newdata = test)
table(correct=correct_classes, predict=lda.pred$class)
## predict
## correct low med_low med_high high
## low 11 6 0 0
## med_low 10 16 3 0
## med_high 1 5 16 2
## high 0 0 0 32
library(MASS)
data("Boston")
boston_scaled<-scale(Boston)
boston_scaled<-as.data.frame(boston_scaled)
dist_eu<-dist(boston_scaled)
summary(dist_eu)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.1343 3.4625 4.8241 4.9111 6.1863 14.3970
km<-kmeans(boston_scaled, centers=3)
pairs(boston_scaled, col=km$cluster)
set.seed(11)
k_max<-20
twcss <- sapply(1:k_max, function(k){kmeans(boston_scaled, k)$tot.withinss})
qqplot(x=1:k_max, y = twcss, geom="Line")
## Warning in plot.window(...): "geom" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "geom" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "geom" is not
## a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "geom" is not
## a graphical parameter
## Warning in box(...): "geom" is not a graphical parameter
## Warning in title(...): "geom" is not a graphical parameter
km<-kmeans(boston_scaled, centers=6)
pairs(boston_scaled, col = km$cluster)
library(GGally)
## Registered S3 method overwritten by 'GGally':
## method from
## +.gg ggplot2
##
## Attaching package: 'GGally'
## The following object is masked from 'package:dplyr':
##
## nasa
ggpairs(boston_scaled)
model_predictors<- dplyr::select(train, -crime)
dim(model_predictors)
## [1] 404 13
dim(lda.fit$scaling)
## [1] 13 3
matrix_product <- as.matrix(model_predictors)%*%lda.fit$scaling
matrix_product <- as.data.frame(matrix_product)
library(plotly)
##
## Attaching package: 'plotly'
## The following object is masked from 'package:MASS':
##
## select
## The following object is masked from 'package:ggplot2':
##
## last_plot
## The following object is masked from 'package:stats':
##
## filter
## The following object is masked from 'package:graphics':
##
## layout
plot_ly(x=matrix_product$LD1, y=matrix_product$LD2, z=matrix_product$LD3, type="scatter3d", mode="markers")
plot_ly(x=matrix_product$LD1, y=matrix_product$LD2, z=matrix_product$LD3, type="scatter3d", mode="markers", color = train$crime)
km <-kmeans(boston_scaled, centers = 6)
plot_ly(x = matrix_product$LD1, y = matrix_product$LD2, z = matrix_product$LD3, type= 'scatter3d', mode='markers', color = km$cluster[ind])
human <- read.table("http://s3.amazonaws.com/assets.datacamp.com/production/course_2218/datasets/human2.txt", sep=",", header=T)
library(MASS);library(tidyr);library(dplyr);library(GGally);library(ggplot2);library(dplyr);library(tidyr);library(corrplot)
str(human);dim(human);summary(human)
## 'data.frame': 155 obs. of 8 variables:
## $ Edu2.FM : num 1.007 0.997 0.983 0.989 0.969 ...
## $ Labo.FM : num 0.891 0.819 0.825 0.884 0.829 ...
## $ Edu.Exp : num 17.5 20.2 15.8 18.7 17.9 16.5 18.6 16.5 15.9 19.2 ...
## $ Life.Exp : num 81.6 82.4 83 80.2 81.6 80.9 80.9 79.1 82 81.8 ...
## $ GNI : int 64992 42261 56431 44025 45435 43919 39568 52947 42155 32689 ...
## $ Mat.Mor : int 4 6 6 5 6 7 9 28 11 8 ...
## $ Ado.Birth: num 7.8 12.1 1.9 5.1 6.2 3.8 8.2 31 14.5 25.3 ...
## $ Parli.F : num 39.6 30.5 28.5 38 36.9 36.9 19.9 19.4 28.2 31.4 ...
## [1] 155 8
## Edu2.FM Labo.FM Edu.Exp Life.Exp
## Min. :0.1717 Min. :0.1857 Min. : 5.40 Min. :49.00
## 1st Qu.:0.7264 1st Qu.:0.5984 1st Qu.:11.25 1st Qu.:66.30
## Median :0.9375 Median :0.7535 Median :13.50 Median :74.20
## Mean :0.8529 Mean :0.7074 Mean :13.18 Mean :71.65
## 3rd Qu.:0.9968 3rd Qu.:0.8535 3rd Qu.:15.20 3rd Qu.:77.25
## Max. :1.4967 Max. :1.0380 Max. :20.20 Max. :83.50
## GNI Mat.Mor Ado.Birth Parli.F
## Min. : 581 Min. : 1.0 Min. : 0.60 Min. : 0.00
## 1st Qu.: 4198 1st Qu.: 11.5 1st Qu.: 12.65 1st Qu.:12.40
## Median : 12040 Median : 49.0 Median : 33.60 Median :19.30
## Mean : 17628 Mean : 149.1 Mean : 47.16 Mean :20.91
## 3rd Qu.: 24512 3rd Qu.: 190.0 3rd Qu.: 71.95 3rd Qu.:27.95
## Max. :123124 Max. :1100.0 Max. :204.80 Max. :57.50
155 oberservations and 8 variables
To see structure and relationships
cor(human)
## Edu2.FM Labo.FM Edu.Exp Life.Exp GNI
## Edu2.FM 1.000000000 0.009564039 0.59325156 0.5760299 0.43030485
## Labo.FM 0.009564039 1.000000000 0.04732183 -0.1400125 -0.02173971
## Edu.Exp 0.593251562 0.047321827 1.00000000 0.7894392 0.62433940
## Life.Exp 0.576029853 -0.140012504 0.78943917 1.0000000 0.62666411
## GNI 0.430304846 -0.021739705 0.62433940 0.6266641 1.00000000
## Mat.Mor -0.660931770 0.240461075 -0.73570257 -0.8571684 -0.49516234
## Ado.Birth -0.529418415 0.120158862 -0.70356489 -0.7291774 -0.55656208
## Parli.F 0.078635285 0.250232608 0.20608156 0.1700863 0.08920818
## Mat.Mor Ado.Birth Parli.F
## Edu2.FM -0.6609318 -0.5294184 0.07863528
## Labo.FM 0.2404611 0.1201589 0.25023261
## Edu.Exp -0.7357026 -0.7035649 0.20608156
## Life.Exp -0.8571684 -0.7291774 0.17008631
## GNI -0.4951623 -0.5565621 0.08920818
## Mat.Mor 1.0000000 0.7586615 -0.08944000
## Ado.Birth 0.7586615 1.0000000 -0.07087810
## Parli.F -0.0894400 -0.0708781 1.00000000
ggpairs(human)
cor(human)%>%corrplot(method="circle", type="upper", title="Graph 1. Correlations between variables in the human data")
pca_human <- prcomp(human)
pca_human
## Standard deviations (1, .., p=8):
## [1] 1.854416e+04 1.855219e+02 2.518701e+01 1.145441e+01 3.766241e+00
## [6] 1.565912e+00 1.912052e-01 1.591112e-01
##
## Rotation (n x k) = (8 x 8):
## PC1 PC2 PC3 PC4
## Edu2.FM -5.607472e-06 0.0006713951 -3.412027e-05 -2.736326e-04
## Labo.FM 2.331945e-07 -0.0002819357 5.302884e-04 -4.692578e-03
## Edu.Exp -9.562910e-05 0.0075529759 1.427664e-02 -3.313505e-02
## Life.Exp -2.815823e-04 0.0283150248 1.294971e-02 -6.752684e-02
## GNI -9.999832e-01 -0.0057723054 -5.156742e-04 4.932889e-05
## Mat.Mor 5.655734e-03 -0.9916320120 1.260302e-01 -6.100534e-03
## Ado.Birth 1.233961e-03 -0.1255502723 -9.918113e-01 5.301595e-03
## Parli.F -5.526460e-05 0.0032317269 -7.398331e-03 -9.971232e-01
## PC5 PC6 PC7 PC8
## Edu2.FM -0.0022935252 2.180183e-02 6.998623e-01 7.139410e-01
## Labo.FM 0.0022190154 3.264423e-02 7.132267e-01 -7.001533e-01
## Edu.Exp 0.1431180282 9.882477e-01 -3.826887e-02 7.776451e-03
## Life.Exp 0.9865644425 -1.453515e-01 5.380452e-03 2.281723e-03
## GNI -0.0001135863 -2.711698e-05 -8.075191e-07 -1.176762e-06
## Mat.Mor 0.0266373214 1.695203e-03 1.355518e-04 8.371934e-04
## Ado.Birth 0.0188618600 1.273198e-02 -8.641234e-05 -1.707885e-04
## Parli.F -0.0716401914 -2.309896e-02 -2.642548e-03 2.680113e-03
biplot(pca_human, choices = 1:2, cex = c(0.8, 1), col= c("grey40", "deeppink2"), main = "Grpagh 2. PCA with human data, unscaled")
## Warning in arrows(0, 0, y[, 1L] * 0.8, y[, 2L] * 0.8, col = col[2L], length
## = arrow.len): zero-length arrow is of indeterminate angle and so skipped
## Warning in arrows(0, 0, y[, 1L] * 0.8, y[, 2L] * 0.8, col = col[2L], length
## = arrow.len): zero-length arrow is of indeterminate angle and so skipped
## Warning in arrows(0, 0, y[, 1L] * 0.8, y[, 2L] * 0.8, col = col[2L], length
## = arrow.len): zero-length arrow is of indeterminate angle and so skipped
## Warning in arrows(0, 0, y[, 1L] * 0.8, y[, 2L] * 0.8, col = col[2L], length
## = arrow.len): zero-length arrow is of indeterminate angle and so skipped
## Warning in arrows(0, 0, y[, 1L] * 0.8, y[, 2L] * 0.8, col = col[2L], length
## = arrow.len): zero-length arrow is of indeterminate angle and so skipped
human.std <- scale(human)
pca_human2 <- prcomp(human.std)
biplot(pca_human2, choices= 1:2, cex = c(0.8, 1), col= c("grey40", "deeppink2"), main="Graph 3. PCA with human data, scaled")
s <- summary(pca_human2)
pca_pr <- round(100*s$importance[2, ], digits = 1)
pc_lab <- paste0(names(pca_pr), " (", pca_pr, "%)")
biplot(pca_human2, cex = c(0.8, 1), col = c("grey40", "deeppink2"), xlab = pc_lab[1], ylab = pc_lab[2], main="Graph 4. PCA with human data, scaled, %")
library(FactoMineR)
es.pca = PCA(human, scale.unit= TRUE, ncp=5)
library(ggplot2)
library(dplyr)
library(tidyr)
library(FactoMineR)
data(tea)
my_tea <- data.frame(tea)
dim(my_tea)
## [1] 300 36
str(my_tea)
## 'data.frame': 300 obs. of 36 variables:
## $ breakfast : Factor w/ 2 levels "breakfast","Not.breakfast": 1 1 2 2 1 2 1 2 1 1 ...
## $ tea.time : Factor w/ 2 levels "Not.tea time",..: 1 1 2 1 1 1 2 2 2 1 ...
## $ evening : Factor w/ 2 levels "evening","Not.evening": 2 2 1 2 1 2 2 1 2 1 ...
## $ lunch : Factor w/ 2 levels "lunch","Not.lunch": 2 2 2 2 2 2 2 2 2 2 ...
## $ dinner : Factor w/ 2 levels "dinner","Not.dinner": 2 2 1 1 2 1 2 2 2 2 ...
## $ always : Factor w/ 2 levels "always","Not.always": 2 2 2 2 1 2 2 2 2 2 ...
## $ home : Factor w/ 2 levels "home","Not.home": 1 1 1 1 1 1 1 1 1 1 ...
## $ work : Factor w/ 2 levels "Not.work","work": 1 1 2 1 1 1 1 1 1 1 ...
## $ tearoom : Factor w/ 2 levels "Not.tearoom",..: 1 1 1 1 1 1 1 1 1 2 ...
## $ friends : Factor w/ 2 levels "friends","Not.friends": 2 2 1 2 2 2 1 2 2 2 ...
## $ resto : Factor w/ 2 levels "Not.resto","resto": 1 1 2 1 1 1 1 1 1 1 ...
## $ pub : Factor w/ 2 levels "Not.pub","pub": 1 1 1 1 1 1 1 1 1 1 ...
## $ Tea : Factor w/ 3 levels "black","Earl Grey",..: 1 1 2 2 2 2 2 1 2 1 ...
## $ How : Factor w/ 4 levels "alone","lemon",..: 1 3 1 1 1 1 1 3 3 1 ...
## $ sugar : Factor w/ 2 levels "No.sugar","sugar": 2 1 1 2 1 1 1 1 1 1 ...
## $ how : Factor w/ 3 levels "tea bag","tea bag+unpackaged",..: 1 1 1 1 1 1 1 1 2 2 ...
## $ where : Factor w/ 3 levels "chain store",..: 1 1 1 1 1 1 1 1 2 2 ...
## $ price : Factor w/ 6 levels "p_branded","p_cheap",..: 4 6 6 6 6 3 6 6 5 5 ...
## $ age : int 39 45 47 23 48 21 37 36 40 37 ...
## $ sex : Factor w/ 2 levels "F","M": 2 1 1 2 2 2 2 1 2 2 ...
## $ SPC : Factor w/ 7 levels "employee","middle",..: 2 2 4 6 1 6 5 2 5 5 ...
## $ Sport : Factor w/ 2 levels "Not.sportsman",..: 2 2 2 1 2 2 2 2 2 1 ...
## $ age_Q : Factor w/ 5 levels "15-24","25-34",..: 3 4 4 1 4 1 3 3 3 3 ...
## $ frequency : Factor w/ 4 levels "1/day","1 to 2/week",..: 1 1 3 1 3 1 4 2 3 3 ...
## $ escape.exoticism: Factor w/ 2 levels "escape-exoticism",..: 2 1 2 1 1 2 2 2 2 2 ...
## $ spirituality : Factor w/ 2 levels "Not.spirituality",..: 1 1 1 2 2 1 1 1 1 1 ...
## $ healthy : Factor w/ 2 levels "healthy","Not.healthy": 1 1 1 1 2 1 1 1 2 1 ...
## $ diuretic : Factor w/ 2 levels "diuretic","Not.diuretic": 2 1 1 2 1 2 2 2 2 1 ...
## $ friendliness : Factor w/ 2 levels "friendliness",..: 2 2 1 2 1 2 2 1 2 1 ...
## $ iron.absorption : Factor w/ 2 levels "iron absorption",..: 2 2 2 2 2 2 2 2 2 2 ...
## $ feminine : Factor w/ 2 levels "feminine","Not.feminine": 2 2 2 2 2 2 2 1 2 2 ...
## $ sophisticated : Factor w/ 2 levels "Not.sophisticated",..: 1 1 1 2 1 1 1 2 2 1 ...
## $ slimming : Factor w/ 2 levels "No.slimming",..: 1 1 1 1 1 1 1 1 1 1 ...
## $ exciting : Factor w/ 2 levels "exciting","No.exciting": 2 1 2 2 2 2 2 2 2 2 ...
## $ relaxing : Factor w/ 2 levels "No.relaxing",..: 1 1 2 2 2 2 2 2 2 2 ...
## $ effect.on.health: Factor w/ 2 levels "effect on health",..: 2 2 2 2 2 2 2 2 2 2 ...
summary(my_tea)
## breakfast tea.time evening lunch
## breakfast :144 Not.tea time:131 evening :103 lunch : 44
## Not.breakfast:156 tea time :169 Not.evening:197 Not.lunch:256
##
##
##
##
##
## dinner always home work
## dinner : 21 always :103 home :291 Not.work:213
## Not.dinner:279 Not.always:197 Not.home: 9 work : 87
##
##
##
##
##
## tearoom friends resto pub
## Not.tearoom:242 friends :196 Not.resto:221 Not.pub:237
## tearoom : 58 Not.friends:104 resto : 79 pub : 63
##
##
##
##
##
## Tea How sugar how
## black : 74 alone:195 No.sugar:155 tea bag :170
## Earl Grey:193 lemon: 33 sugar :145 tea bag+unpackaged: 94
## green : 33 milk : 63 unpackaged : 36
## other: 9
##
##
##
## where price age sex
## chain store :192 p_branded : 95 Min. :15.00 F:178
## chain store+tea shop: 78 p_cheap : 7 1st Qu.:23.00 M:122
## tea shop : 30 p_private label: 21 Median :32.00
## p_unknown : 12 Mean :37.05
## p_upscale : 53 3rd Qu.:48.00
## p_variable :112 Max. :90.00
##
## SPC Sport age_Q frequency
## employee :59 Not.sportsman:121 15-24:92 1/day : 95
## middle :40 sportsman :179 25-34:69 1 to 2/week: 44
## non-worker :64 35-44:40 +2/day :127
## other worker:20 45-59:61 3 to 6/week: 34
## senior :35 +60 :38
## student :70
## workman :12
## escape.exoticism spirituality healthy
## escape-exoticism :142 Not.spirituality:206 healthy :210
## Not.escape-exoticism:158 spirituality : 94 Not.healthy: 90
##
##
##
##
##
## diuretic friendliness iron.absorption
## diuretic :174 friendliness :242 iron absorption : 31
## Not.diuretic:126 Not.friendliness: 58 Not.iron absorption:269
##
##
##
##
##
## feminine sophisticated slimming
## feminine :129 Not.sophisticated: 85 No.slimming:255
## Not.feminine:171 sophisticated :215 slimming : 45
##
##
##
##
##
## exciting relaxing effect.on.health
## exciting :116 No.relaxing:113 effect on health : 66
## No.exciting:184 relaxing :187 No.effect on health:234
##
##
##
##
##
gather(my_tea) %>% ggplot(aes(value)) + facet_wrap("key", scales = "free") + geom_bar() + theme(axis.text.x = element_text(angle = 45, hjust = 1, size = 8))
## Warning: attributes are not identical across measure variables;
## they will be dropped
keep_columns <- c("breakfast", "sex", "price", "healthy", "spirituality", "friends", "Tea", "tearoom")
my_tea1 <- dplyr::select(my_tea, one_of(keep_columns))
gather(my_tea1) %>% ggplot(aes(value)) + facet_wrap("key", scales = "free") + geom_bar() + theme(axis.text.x = element_text(angle = 45, hjust = 1, size = 8))
## Warning: attributes are not identical across measure variables;
## they will be dropped
mca <- MCA(my_tea1, graph = FALSE)
summary(mca)
##
## Call:
## MCA(X = my_tea1, graph = FALSE)
##
##
## Eigenvalues
## Dim.1 Dim.2 Dim.3 Dim.4 Dim.5 Dim.6
## Variance 0.193 0.165 0.150 0.147 0.139 0.128
## % of var. 11.885 10.164 9.218 9.031 8.544 7.859
## Cumulative % of var. 11.885 22.050 31.268 40.299 48.843 56.702
## Dim.7 Dim.8 Dim.9 Dim.10 Dim.11 Dim.12
## Variance 0.121 0.115 0.110 0.102 0.096 0.087
## % of var. 7.429 7.054 6.784 6.297 5.928 5.371
## Cumulative % of var. 64.131 71.185 77.970 84.267 90.195 95.565
## Dim.13
## Variance 0.072
## % of var. 4.435
## Cumulative % of var. 100.000
##
## Individuals (the 10 first)
## Dim.1 ctr cos2 Dim.2 ctr cos2 Dim.3
## 1 | 0.767 1.015 0.144 | 0.347 0.242 0.029 | 0.713
## 2 | 0.261 0.117 0.057 | 0.285 0.164 0.069 | -0.327
## 3 | -0.400 0.276 0.233 | -0.236 0.113 0.081 | 0.005
## 4 | 0.057 0.006 0.003 | -0.199 0.080 0.034 | 0.209
## 5 | -0.005 0.000 0.000 | 0.067 0.009 0.003 | 0.484
## 6 | 0.641 0.709 0.171 | -0.080 0.013 0.003 | 0.234
## 7 | -0.151 0.039 0.028 | 0.019 0.001 0.000 | 0.020
## 8 | 0.262 0.118 0.059 | 0.161 0.052 0.022 | -0.149
## 9 | 0.533 0.491 0.180 | 0.465 0.437 0.137 | 0.375
## 10 | 0.503 0.436 0.118 | 1.130 2.577 0.594 | -0.477
## ctr cos2
## 1 1.132 0.125 |
## 2 0.238 0.090 |
## 3 0.000 0.000 |
## 4 0.097 0.037 |
## 5 0.521 0.164 |
## 6 0.122 0.023 |
## 7 0.001 0.001 |
## 8 0.050 0.019 |
## 9 0.313 0.089 |
## 10 0.505 0.106 |
##
## Categories (the 10 first)
## Dim.1 ctr cos2 v.test Dim.2 ctr
## breakfast | -0.002 0.000 0.000 -0.030 | 0.211 1.613
## Not.breakfast | 0.002 0.000 0.000 0.030 | -0.194 1.489
## F | -0.357 4.897 0.186 -7.458 | -0.173 1.347
## M | 0.521 7.144 0.186 7.458 | 0.253 1.965
## p_branded | 0.289 1.708 0.039 3.398 | -0.960 22.087
## p_cheap | 0.678 0.694 0.011 1.812 | 0.329 0.191
## p_private label | 0.723 2.369 0.039 3.431 | 0.274 0.398
## p_unknown | 0.192 0.095 0.002 0.677 | -0.040 0.005
## p_upscale | 0.563 3.621 0.068 4.507 | 1.183 18.709
## p_variable | -0.710 12.168 0.300 -9.471 | 0.187 0.986
## cos2 v.test Dim.3 ctr cos2 v.test
## breakfast 0.041 3.500 | -0.286 3.276 0.075 -4.751 |
## Not.breakfast 0.041 -3.500 | 0.264 3.024 0.075 4.751 |
## F 0.044 -3.618 | -0.242 2.905 0.086 -5.059 |
## M 0.044 3.618 | 0.353 4.238 0.086 5.059 |
## p_branded 0.427 -11.300 | -0.341 3.069 0.054 -4.012 |
## p_cheap 0.003 0.880 | 1.264 3.110 0.038 3.378 |
## p_private label 0.006 1.301 | 0.150 0.131 0.002 0.709 |
## p_unknown 0.000 -0.142 | 2.649 23.430 0.292 9.351 |
## p_upscale 0.300 9.475 | -0.267 1.049 0.015 -2.137 |
## p_variable 0.021 2.493 | 0.024 0.019 0.000 0.326 |
##
## Categorical variables (eta2)
## Dim.1 Dim.2 Dim.3
## breakfast | 0.000 0.041 0.075 |
## sex | 0.186 0.044 0.086 |
## price | 0.319 0.560 0.369 |
## healthy | 0.010 0.044 0.413 |
## spirituality | 0.082 0.019 0.000 |
## friends | 0.412 0.000 0.000 |
## Tea | 0.272 0.340 0.162 |
## tearoom | 0.264 0.273 0.092 |
plot(mca, invisible=c("ind"), habillage="quali")
res.mca=MCA(my_tea1)